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For  several  years,  SRI  hasaconducted/S  program  under  sponsorship  of  the 

Defense  Nuclear  Agency  (DNA)  to  measure4the  dynamic  mechanical  response  of  in 
situ  geologic  materials  near  and  above  the  elastic  limit.  A major  achievement- 
in  this  program  '.has"  been  development  of  the  LASS  technique  (Lagrange  Analysis 
for  Stresses  and  Strains),  in  which  a mathematically  complete  Lagrange  flow 
analysis  is  used  to  calculate  dynamic  stress-strain  trajectories  from  field 
measurements  of  explosive-induced  ground  motion.  , , 
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In  the  present  work  three  tasks  were  completed:  (1)  two  fields  shots 
were  fired  to  provide  data  for  proof  testing  the  LASS  technique?  (2)N  the  field 
test  data  were  reduced  using  LASS  I and^-f^^LASS  was  extended  and  its  accuracy 
and  sensitivity  to  typical  errors  in  field  measurements  were  evaluated. 

The  two  field  tests  were  116-kg  spherical  shots  fired  in  a wet  clay  layer 
about  3.6  m below  the  surface  at  the  Pre-Dice  Throw  11  site.  Excellent  stress 
and  particle  velocity  records  were  obtained  with  peak  stress  between  0.2  and 
1.7  GPa.  , Previously  developed  instrumentation  methods  were  improved  in  the 
course  of)  these  tests. 

^Plausible  stress-strain  trajectories  were  calculated  from  the  field 
records  using  LASS.  These  indicate  nearly  reversible  radial  stress-volume 
strain  behavior  unt il jrelease  to  low  stresses,  possible  bulking  at  low  stresses 
on  release,  and  about  one-third  the  compressibility  observed  in  laboratory 
uniaxial  strain  tests.  The  calculated  principal  stress  difference  amplitude 
|was  small,  less  than  the  experimental  error  associated  with  the  calculation. 
This  result  was  expected  because  of  the  high  stress  level  of  the  tests  and  the 
low  strength  of  the  site  material.  "f— 

LASS  analysis  capabilities  were  extended  by  the  addition  of  variable 
window-width,  adjustable  low-pass  digital  smoothing  procedures  for  the  field 
records,  the  derivation  of  statistical  error  propagation  equations  for  the 
stress  and  strain  calculations,  and  the  formulation  and  reduction  to  code  of  a 
new  and  more  flexible  form  of  LASS  called  the  global  Lagrange  analysis. 

Analysis  validation  calculations  performed  on  constructed  problems  for  which 
analytic  or  numerical  solutions  are  available  showed  that,  with  a realistic 
number  of  gage  stations  (five)  at  which  perfect  stress  and  particle  velocity 
records  are  obtained:  (1)  stress  difference  calculations  are  accurate  to 
|within  10%,  if  the  peak  stresses  at  the  measurement  stations  are  comparable  to 
the  stress  difference  and  (2)  strain  calculations  are  accurate  to  within  5%. 
Sensitivity  studies  showed  that  piecewise  linear  and  global  forms  of  LASS 
remain  stable  when  small  amounts  of  random  amplitude  or  position  noise  are 
added  to  the  gage  records. 

Recommendations  for  four  possible  future  research  activities  related  to 
the  present  results  are  offered:  (1)  experiments  to  improve  the  reliability 
lof  in  situ  stress  measurements  in  soils  at  low  stresses  where  soil  strength 
effects  presumably  preclude  the  use  of  the  techniques  employed  in  the  present 
effort,  (2)  the  addition  of  LASS  capabilities  to  CIST,  the  Air  Force  Weapons 
(Laboratory  standardized  cylindrical  in  situ  test,  (3)  sensitivity  studies  and 
error  analyses  using  LASS  to  support  and  optimize  the  design  of  proposed  buried 
in  situ  shots  to  determine  material  properties,  and  (4)  theoretical  studies  to 
evaluate  the  sensitivity  of  LASS  to  deviations  from  the  assumed  flow  symmetry. 
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SUMMARY 


BACKGROUND 


The  overall  objectives  of  this  multiyear  program  have  been  to 
* 

develop  the  LASS  method  for  determining  in  situ  mechanical  response 
of  geologic  materials  and  to  apply  it  to  field  data  generated  at  sites 
of  interest  to  DNA.  Previous  accomplishments  in  the  program  (see  Appendix  B) 
include  development  of  the  ytterbium  stress  gage  and  the  mutual  inductance 
particle  velocity  gage  for  rock  and  soil  measurements  at  high  amplitudes 
and  loading  rates,  application  of  these  gages  for  close-in  measurements 
on  events  of  concern  to  DNA,  and  development  of  LASS.  The  objectives 
of  the  present  work  were  to  modify  experimental  techniques  and  the  Lagrange 
analysis  (LASS)  for  increased  accuracy  and  to  proof  test  LASS  by  applying 
it  to  data  generated  in  field  tests  at  a site  selected  jointly  by  DNA  and 
SRI. 


The  LASS  method  consists  of  applying  a mathematical  Lagrange  flow 
analysis  to  in  situ  measurements  of  radial  stress  and  radial  particle 
velocity  made  during  explosive-induced  spherical  or  cylindrical  site 
deformation.  Results  are  the  stress  and  strain  trajectories  defining  in 
situ  response  to  the  applied  large-amplitude  dynamic  deformation.  A 
key  point  is  that  for  flows  satisfying  the  spherical  or  cylindrical 
symmetry  condition,  the  analysis  in  conjunction  with  the  field  measurements 
is  mathematically  complete.  Therefore  (1)  the  LASS  results  are  unique 
in  the  sense  that  no  assumptions  concerning  the  constitutive  relations 
of  the  test  material  are  used  in  the  analysis,  and  (2)  all  independent 
stresses^  and  strains  are  calculable  for  particles  within  the  instrumented 
region  of  the  flow. 
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Lagrange  Analysis  for  Stress  and  Strains. 

^Axial  stress  is  not  obtained  for  cylindrical  flow  because  it  does  not 
appear  in  the  conservation  of  radial  momentum  law  for  that  geometry. 
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An  in  situ  dynamic  loading  test  known  as  CIST  has  been  developed 
by  the  Air  Force  Weapons  Laboratory,  Because  it  employs  a cylindrical 
loading  source,  CIST  can  simultaneously  test  several  strata  within  a 
layered  site  and  hence  is  cost-effective  for  such  a site.  However, 
because  CIST  data  are  analyzed  by  an  iterative  procedure  in  which  a model 
is  assumed  and  the  free  parameters  adjusted  until  satisfactory  agreement 
with  displacement-related  field  measurements  is  obtained,  the  resulting 
material  properties  are  not  unique  and  may  be  incomplete  because  of 
inadequacies  in  the  assumed  model.  For  site  regions  in  which  the  CIST 
source  generates  cylindrical  flow,  we  suggest  combining  the  LASS  analysis 
with  CIST  field  techniques,  augmented  by  stress  measurements,  to  obtain 
a site  characterization  method  of  high  accuracy  and  cost-effectiveness. 

PROGRESS 

Progress  in  the  LASS  program  during  the  past  year  is  summarized 
below  under  the  following  three  headings: 

• Field  tests 

• LASS  analyses  of  the  field  tests 

• Analysis  extensions  and  validation. 

Field  tests 

The  field  tests  consisted  of  two  replicate  spherical  shots  performed 
in  the  Queen  15  Area  of  the  White  Sands  Missile  Range  (WSMR) . The  shot 
sites  were  within  about  200  m of  both  Pre-Dice  Throw  II  events  and  CIST  15 
(see  Figure  1 in  text) . The  test  layer  was  a relatively  uniform  wet  clay 
stratum  between  the  water  table  at  a depth  of  about  2.1  m and  a more  sandy 
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Cylindrical  In  Situ  Test. 
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clay  layer  at  about  5.2  m.  The  HE  sources  were  116-kg  cast  TNT  spheres. 
Instrumentation  layout  and  dimensions  are  shown  in  Figures  2 through  4. 

Two  significant  instrumentation  improvements  were  incorporated  in 
these  tests.  First,  the  mutual  inductance  particle  velocity  gage  power 
supplies  were  modified  to  permit  remote  operation.  By  locating  the 
power  supplies  near  the  shot,  we  significantly  reduced  noise  generation 
and  pickup  problems,  power  losses  in  the  current  cable,  and  cabling 
costs.  Second,  electrical  cross  talk  between  gages  was  reduced  by 
modifying  locations  and  orientations  of  gages  and  cables  and  by  avoiding 
ground  loops. 

The  shots  yielded  a 100%  data  return  consisting  of  exceptionally 
high  quality  stress  and  particle  velocity  records.  Peak  stresses  were 
between  0.25  and  1.7  GPa,  and  peak  particle  velocities  were  between 
50  and  175  ms  . The  reduced  records  are  shown  in  Figures  8 through  15. 

LASS  Analysis  of  Field  Tests 

Stress  and  strain  trajectories  for  the  test  layer  were  calculated 
from  the  field  data  using  LASS.  The  calculated  trajectories,  shown  in 
Figures  27  through  29  and  37  through  39,  are  estimated  to  be  accurate  to 
about  20%  of  their  peak  values.  The  agreement  between  the  two  tests, 
as  shown  in  Figure  41,  is  about  20%  based  on  the  average  compressibility 
calculated  at  1.05  m,  the  middle  of  the  measurement  range.  This  is  well 
within  both  experimental  error  and  the  range  of  variation  expected  for 
this  geological  site.  The  material  response  characteristics  indicated 
by  the  calculated  dynamic  stress-strain  trajectories,  including  no 
hysteresis  in  the  load-release  cycle  and  no  irreversible  crushing,  are 
quite  reasonable  for  wet  clay  with  virtually  no  unfilled  voids.  Because 
of  the  high  stress  level  of  the  field  tests  compared  with  the  stress 
difference  $ supportable  by  the  test  layer,  the  errors  accruing  in  the 
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$ calculation  are  greater  than  $ everywhere  but  during  loading.  Therefore, 
to  the  precision  of  these  tests,  may  be  taken  to  be  zero.  Comparison 

of  dynamic  loading  states  measured  in  the  proof  tests  with  laboratory 
uniaxial  strain  test  data  obtained  by  the  U.S.  Army  Waterways  Experiment 
Station  for  similar  material  from  the  Pre-Dice  Throw  II  site  indicates 
qualitatively  similar  behavior  but  about  two-thirds  less  compressibility 
in  the  LASS  tests  (Figure  41) . Thus  the  wet  clay  layer  appears  to  be 
sensitive  to  the  differences  between  laboratory  and  in  situ  test  conditions. 

The  proof  tests  are  judged  to  have  been  very  successful  in  that 
they  validated  the  LASS  technique  and  also  provided  material  properties 
for  the  Pre-Dice  Throw  II  site.  We  conclude  that  LASS  is  ready  for 
application  at  appropriate  sites  of  interest  to  DNA, 


Analysis  Extensions  and  Validation 

We  made  three  major  extensions  to  LASS  capabilities  in  the  present 
work.  First,  digital  data  smoothing  techniques  based  on  modified  low-pass 
filtering  concepts  were  added  to  maximize  uniformity  and  reproducibility 
in  data  reduction.  Second,  an  alternative  Lagrange  analysis,  termed  the 
global  form,  was  formulated.  This  form  makes  more  nearly  optimal  use  of 
the  field  data,  is  more  convenient  for  error  propagation  analyses,  provides 
built-in  data  smoothing,  and  permits  enforcement  of  physical  constraints, 
if  desired.  Third,  we  derived  expressions  for  the  propagated  errors  in 
calculated  stresses  and  strains  for  each  form  of  the  analysis.  We  used 
standard  statistical  methods  with  first-order  mathematical  approximations 
and,  when  applicable,  simplifying  assumptions.  For  the  global  Lagrange 
analysis,  the  derived  expressions  were  reduced  to  code  form. 

To  validate  the  analysis  we  performed  calculations  with  either  the 
piecewise  linear  or  global  forms  of  LASS  on  an  elastic  test  problem  for 
which  we  have  analytic  solutions,  on  an  inelastic  test  problem  that  we 


4 


created  last  year,  and  on  LASS  1 data.  As  a result  we  made  two  technique 
modifications  to  the  piecewise  linear  analysis,  resulting  in  significant 
improvements  in  the  accuracy  of  calculated  strains  and  stress  difference. 

For  the  elastic  problem  with  five  gage  stations,  we  achieved  10%  or  better 
accuracies  for  all  quantities  including  $ in  the  center  of  the  measurement 
interval  (Figure  44a) . For  the  more  complex  inelastic  problem,  comparable 
accuracies  were  obtained  for  the  strains  (Figure  45) , but  the  calculated 
histories  were  noisy,  primarily  because  $ is  small  compared  with  the 
radial  stresses  in  this  flow  region. 

The  global  analysis  was  applied  to  the  elastic  problem  and  demonstrated 
to  work  properly.  It  was  also  applied  to  the  same  problem  after  random 
noise  and  gage  location  errors  were  added.  The  solutions  were  not 
significantly  degraded  and  were  consistent  with  the  error  propagation 
analysis  predictions.  We  also  demonstrated  that  reducing  the  number  of 
gage  stations  from  five  to  three  approximately  doubles  the  analysis 
errors.  (These  results  are  shown  in  Figure  46  cases  3 and  4,  last  two 
plots.)  A crude  global  analysis  of  the  first  proof  test  provided 
approximate  stress-strain  paths  and  more  importantly  an  estimate  of  the 
associated  errors  resulting  from  measurement  uncertainties  (see  Figure 
47)  . 

We  conclude  from  the  validation  studies  that  both  forms  of  the 
Lagrange  analysis  are  working  properly,  that  we  can  estimate  propagated 
errors  using  the  global  analysis,  that  strains  can  be  accurately  calculated 
from  a realistic  number  of  high  quality  gage  records,  and  that  the  stress 
difference  can  be  accurately  calculated  for  flows  in  which  it  is  comparable 
in  amplitude  to  the  measured  stresses. 

FUTURE  WORK 

For  future  work  we  recommend  (1)  an  experimental  program  to  improve 
in  situ  dynamic  stress  measurement  capabilities  in  the  important  range 
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between  10  and  10  GPa  where  properties  are  needed  and  are 
currently  not  well  understood,  (2)  the  application  of  LASS  analysis 
techniques  to  appropriate  CISTs,  if  reliable  stress  measurements  can  be 
made,  (3)  pretest  LASS  analysis  and  error  propagation  calculations  to 
support  the  design  of  in  situ  tests  intended  to  provide  material  property 
results  of  a specified  accuracy,  and  (4)  sensitivity  studies  to  evaluate 
the  magnitude  of  errors  caused  by  applying  LASS  methods  to  data  from 
slightly  nonuniform  or  anisotropic  site  layers. 
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1.  INTRODUCTION 


For  several  years  the  Defense  Nuclear  Agency  (DNA)  has  conducted 
a program  at  SRI  to  determine  dynamic  in  situ  constitutive  properties 
of  soils  and  rocks  at  sites  of  interest.  The  program  has  had  three 
major  emphases: 

• Development  of  ground  motion  instrumentation  in  the 
stress  range  from  0.1  to  3.0  GPa. 

• Development  of  a one-dimensional  divergent  flow 
Lagrange  analysis. 

• Field  testing. 

The  program  has  produced  a new  site  characterization  method  that  draws 
on  results  from  each  of  the  above  areas.  In  this  technique  a Lagrange 
analysis  is  applied  to  field  measurements  of  site  material  motion  induced 
by  explosive  loading  to  obtain  dynamic  stress-strain  trajectories  for 
the  in  situ  material.  We  refer  to  the  reduction  of  appropriate  field 
data  by  Lagrange  analysis  as  the  LASS  technique  (Lagrange  analysis  for 
stresses  and  strains) . This  report  describes  new  material  properties 
results  obtained  with  LASS  and  recent  developments  of  MSS  capabilities. 

PROBLEM 

Determination  of  the  vulnerability  of  structures  to  near  surface 
detonations  involves  calculating  the  airblast,  ground  shock,  cratering, 
and  soil-structure  interactions  induced  by  the  event.  However,  for 
partially  contained  detonations  and  for  complicated  site  geologies  and 
materials,  code  predictions  of  these  phenomena  often  do  not  agree  with 
field  measurements.  The  discrepancies  are  usually  attributed  to 
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inadequacies  in  the  material  models  or  the  material  properties  used  to 
represent  the  site  material  in  the  calculations.  Properties  of  the  yield 
surface  are  particularly  uncertain. 

For  HE  events  the  dynamic  mechanical  response  for  in  situ  material, 

along  a variety  of  loading  paths  in  stress  space,  for  stresses  from  under 

-2  -1 

a MPa  to  well  over  a GPa,  and  at  strain  rates  from  less  than  10  s to 
4 -1 

greater  than  10  s must  be  correctly  specified  by  the  material  model  if 
soil  response  is  to  be  satisfactorily  calculated.  Traditionally,  a fairly 
extensive  data  base  for  constructing  these  models  would  contain  information 
on  generic  material  properties  from  previous  studies  on  presumably  similar 
materials,  the  results  of  various  mechanical  and  physical  tests  performed 
in  the  laboratory  on  cores  removed  from  the  test  site,  and  some  on-site 
characterization  data  such  as  the  results  of  seismic  and  geologic  surveys. 
Unfortunately,  since  such  data  are  not  always  obtained  under  conditions 
matching  those  existing  in  the  event  to  be  calculated,  a model  constructed 
from  them  could  possibly  be  seriously  in  error  when  applied  to  the  event. 
Specific  factors  that  are  different  from  the  event  conditions  in  one  or 
more  of  the  data  base  sources  listed  above,  and  that  are  capable  of  signif- 
icantly affecting  specimen  material  response,  include  (1)  in  situ  environ- 
mental conditions  such  as  specimen  moisture  content  and  preexisting  litho- 
static  stresses,  (2)  specimen  crack  and  defect  content  calculated  over 
volumes  scaled  to  the  predominant  wavelengths  expected  in  the  event  being 
calculated,  and  (3)  loading  conditions  including  amplitude,  rate,  and 
stress  path. 

Various  computational  and  experimental  approaches  have  been  pursued, 
(a)  to  determine  the  relative  importance  of  the  present  limitations  in 
the  material  properties  data  base,  (b)  to  provide  "correction  factors" 
that  account  for  phenomenology  not  adequatly  addressed  to  date  (e.g., 
rate  effects,  anisotropy,  etc.)  and  (c)  to  develop  testing  equipment 
and  techniques  to  address  these  shortcomings.  The  approaches  include 
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code  studies  of  parameter  sensitivity;  laboratory  rate  dependence  and 
prestress  tests,  usually  with  small  specimens;  in  situ  quasi-static 
flat  jack  tests;  the  use  of  larger,  faster  laboratory  mechanical  testing 
devices;  and,  recently,  iterative  model  adjustments  based  on  field 
measurements  (CIST).  In  spite  of  such  efforts,  the  constitutive 
relations  for  dynamically  loaded  in  situ  material  remain  a major  uncer- 
tainty in  ground  motion,  cratering,  soil-structure  interaction,  and 
vulnerability  calculations. 

Ideally,  what  is  needed  is  site  material  response  information  obtained 
under  field  event  conditions.  In  principle,  tests  should  be  performed 
with  representative  specimens  of  site  material  comparable  in  size  to 
the  wavelengths  expected  in  the  actual  event,  in  the  in  situ  environment, 
at  the  loading  rates  and  amplitudes  of  interest  in  the  event,  and  along 

jfc 

several  different  stress  paths. 

APPROACH 

The  method  SRI  has  developed  for  determining  site  constitutive 
properties  from  which  to  construct  material  models  is  called  the  LASS 
technique,  after  the  Lagrange  analysis  used  to  reduce  the  data.  In  this 
technique,  a buried  high  explosive  (HE)  charge,  usually  100  kg  or  larger, 
is  used  to  induce  either  spherical  or  cylindrical  flow  over  several  cubic 
meters  or  more  of  the  in  situ  site  material  of  interest.  Radial  stress 
and  radial  particle  velocity  histories  along  several  particle  paths  in 
the  induced  flow  are  measured  with  appropriate  ground  motion  gages, 

l 

♦ 

Since  no  single  experiment  is  capable  of  providing  all  of  this  information 
for  completely  characterizing  a given  site,  it  is  essential  that  the 
various  available  techniques  be  correlated  and  integrated  to  produce  an 
optimal,  cost-effective,  site  characterization.  Field  experiments 
comprise  a key  element  in  this  approach. 
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and  the  field  measurements  are  used  in  a Lagrange  flow  analysis  to 
calculate  the  principal  stress  and  strain  histories  of  particles  at 
various  locations  in  the  flow.  The  resulting  stress  and  strain  histories 
are  an  important  addition  to  the  data  base  for  material  modeling  for 
buried  events  because  they  apply  under  conditions  close  to  those  existing 
in  the  events  to  be  calculated,  they  are  readily  related  to  model 
parameters,  and  they  can  be  used  to  modify  or  validate  previously 
constructed  models. 

The  Lagrange  analysis  has  been  described  in  detail  in  previous 
1-3 

reports  and  is  summarized  in  Appendix  A.  Basically  it  is  a numerical 
analysis  in  which  field  measurements  of  radial  stress  and  radial  particle 
velocity  along  several  particle  paths  are  used  to  estimate  partial 
derivatives  appearing  in  the  conservation  of  mass  and  momentum  laws  for 
one-dimensional  divergent  flows  (spherical  or  cylindrical) . The  conser- 
vation laws  can  then  be  solved  for  the  material  density  and  the  difficult- 
to-measure  tangential  principal  stress  as  functions  of  particle  coordinate 
and  time.  The  principal  engineering  strains  as  functions  of  particle 
position  and  time  are  evaluated  from  their  definitional  equations  and 

the  measurements  of  radial  particle  velocity.  Since  the  analysis  determines 

* 

all  the  principal  stresses  and  strains  as  function  of  Lagrange  coordinate 
and  time,  any  desired  stress-strain  cross-plot  for  a given  particle 
within  the  instrumented  region  of  the  flow  is  available  by  eliminating 
the  time  parameter  at  that  Lagrange  coordinate. 

The  LASS  technique  is  well  suited  to  providing  input  for  the  material 
models  of  interest  to  DNA  because  it  is  applied  to  data  obtained  in  moder- 
ately large  HE  field  tests  and  because  it  uses  a mathematically  complete 
analysis  to  reduce  the  data.  The  advantages  of  obtaining  data  in  a 


*In  cylindrical  flow  the  axial  stress  is  not  determined  since  it  does 
not  appear  in  the  conservation  equations. 
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large-scale  HE  field  test  are  that  the  results  apply  to  in  situ  material, 
and  that  the  specimen  volumes  are  great  enough  to  statistically  average 
the  effects  of  many  in  situ  flaws,  defects,  and  inclusions  that  would 
not  be  present  in  the  smaller  laboratory  specimens.  The  loading  rates 
are  greater  than  those  usually  achieved  in  the  laboratory  and  provide 
an  upper  bound  for  (or  in  some  cases  are  comparable  to)  those  occurring 
in  DNA  problems.  These  advantages  also  apply,  of  course,  to  any  similar 
large  amplitude  dynamic  in  situ  test,  such  as  CIST,  discussed  later  in 
this  section. 

Since  the  Lagrange  analysis  is  mathematically  complete,  it  has 

several  advantages  over  in  situ  material  properties  determination  that 

do  not  have  this  property.  The  results  are  unique  in  the  sense  that  no 

assumptions  concerning  the  constitutive  relations  of  the  test  material, 

and  no  iterative  procedures,  are  used  in  data  reduction.  Also,  all  the 

* 

independent  stresses  and  strains,  including  much-needed  shear  information, 
can  be  calculated  for  particles  within  the  instrumented  region  of  the 
flow.  Furthermore,  LASS  output  is  stress  and  strain  trajectories — a 
convenient  form  for  constructing  material  models;  error  estimates  to 
establish  confidence  limits  for  the  trajectories  can  be  made;  and  the 
technique  is  applicable  in  any  stress  range  or  material  response  regime 
as  long  as  the  flow  is  one-dimensional  and  the  required  measurements 
are  obtained. 

The  limitations  of  LASS  for  site  modeling  are  of  two  types:  field 
data  limitations  and  analysis  limitations.  Field  data  obtained  by  any 
means  are  usually  noisier  and  less  reliable  than  laboratory  data  both 
because  of  the  uncontrolled  environment  in  which  the  measurements  are 
made  and  because  of  the  likelihood  of  undetected  natural  inhomogeneities 


* 

Except  for  axial  stress  in  cylindrical  flows. 
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within  the  test  material.  In  addition,  emplacement  of  the  source  and/or 
instrumentation  may  alter  the  site  properties  both  in  the  immediate 
vicinity  of  the  transducers  and  throughout  the  test  region. 

The  LASS  analysis  limitations  are  that  (1)  it  is  strictly  valid 
only  within  homogeneous  isotropic  site  regions  that  support  one-dimensional 
flows,  and  (2)  the  accuracy  of  the  results  depends  on  both  the  quality 
of  the  field  data  and  the  number  of  stations  at  which  the  data  are 
obtained.  Neither  limitation  is  considered  severe.  The  homogeneity  and 
isotropy  constraints  are  consistent  with  conditions  observed  at  many 
sites  and  assumed  in  many  of  the  models  and  calculations  that  LASS  is 
intended  to  support.  Accuracy  of  any  test  is  limited  by  the  quality  and 
quantity  of  the  measurements;  however,  the  present  work  demonstrates 
that  LASS  produces  satisfactory  precision  with  realistic  numbers  of 
field  quality  records  a-d  is  not  unduly  sensitive  to  the  uncertainties 
included  in  field  measurements. 

The  Air  Force  Weapons  Laboratory  (AFWL)  has  developed  a cylindrical 
4 

in  situ  test,  CIST,  which  has  provided  important  ground  motion  information 
at  a number  of  sites  of  interest  to  DNA.  In  CIST  a relatively  standardized 
cylindrical  HK  source  is  used  to  load  the  site,  and  the  resulting  horizontal, 
and  sometimes  vertical,  acceleration  histories  are  measured  in  identified 
site  layers  at  various  ranges.  Cavity  pressures  are  also  measured  and 
usually  permit  an  idealized  representation  of  the  source  in  that  test. 
Attempts  at  soil  stress  measurement  hove  so  far  not  been  successful 
because  of  experimental  difficulties,  such  as  bridging  aiound  the  gage, 
encountered  at  the  low  stress  levels  (usually  less  than  50  MPa)  examined 
in  a typical  CIST. 

Material  properties  are  calculated  from  CIST  data  by  an  iterative 
process  in  which  a model  is  assumed  and  the  parameters  are  adjusted 
until  the  best  agreement  with  the  measured  particle  velocity  (integrated 


acceleration)  histories  is  obtained.  Although  this  analysis  does  not 
produce  unique  results,  it  is  necessary  because  the  flow  is  not  fully 
specified  by  velocity  histories  alone.  If  stress  histories  were  available, 
much  stronger  constraints  would  be  put  on  the  model.  In  some  cases,  such 
as  for  cylindrical  flows,  the  LASS  analysis  could  be  used.  The  results 
would  provide  a baiis  for  constructing  the  models  before  any  iterations 
and  would  significantly  increase  the  accuracy  and  reliability  of  the 
resulting  constitutive  relations. 

Because  of  its  cylindrical  geometry,  CIST  is  a cost  effective  means 
for  examining  site  properties  as  a function  of  depth.  It  is  also  well 
suited  to  layered  sites.  If  the  layers  are  homogeneous  and  isotropic 
so  the  flow  is  strictly  cylindrical  within  layers  until  communication 
from  the  boundaries  occurs,  the  LASS  analysis  may  be  used  for  data 
reduction  if  stress  measurements  are  obtained.  After  layer  interactions, 
or  for  transversely  isotropic  layers,  CIST  still  provides  information 
for  evaluating  the  parameters  in  appropriate  models;  however,  the 
uniqueness  problem,  already  significant  in  homogeneous  isotropic  noninteracting 
layers,  becomes  severe.  It  probably  would  not  be  feasible  to  use  CIST  to 
examine  or  model  higher  degrees  of  site  anisotropy;  nor  are  such  levels 
of  complexity  now  likely  to  be  retained  in  either  material  models  or 
DNA  ground  motion  calculations.  For  appropriate  sites  we  have  recommended 
that  the  LASS  analysis  be  used  in  conjunction  with  CIST's,  if  the  experi- 
mental problems  in  measuring  low  stresses  in  soils  and  rocks  are  solved, 
to  obtain  a site  characterization  method  of  high  accuracy  and  cost 
effectiveness. 

PREVIOUS  RESULTS 

Previous  work  in  the  SRI  in  situ  properties  program  and  related 

i-3,5-14 

areas  has  been  described  in  a series  of  DNA  reports  summarized 

in  Appendix  B.  Major  milestones  were  the  proposal  of  the  LASS  concept 
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and  the  demonstration  of  its  feasibility,  the  development  of  ytterbium 

stress  gages  and  mutual  inductance  particle  velocity  gages  to  the  status 

2 5-7 

of  working  field  gages,  ’ the  subsequent  use  of  ytterbium  stress 

gages  to  make  ground  motion  measurements  for  DNA  in  a number  of  nuclear 
8,10-14 

and  HE  events,  the  development  of  the  piecewise  linear  LASS 

analysis  (see  Appendix  A)  to  code  form  in  which  it  can  calculate  stress- 

1-3  5-7 

strain  trajectories  for  complex  time-dependent  flows,  ' ’ demonstration 

that  the  accuracy  of  the  numerical  LASS  methods  for  calculating  strains 

when  applied  to  perfect  records  from  5 to  8 gage  stations  is  better 

than  10%  , and  the  first  field  application  of  LASS  in  the  WES  CENSE  II 
1 

series . 

REPORT  ORGANIZATION 

Section  2 describes  the  layout,  instrumentation,  and  data  from  this 
year's  field  tests.  Section  3 contains  the  LASS  analysis  of  the  field 
tests,  a discussion  of  the  results,  and  a comparison  with  other  data. 
Section  4 describes  analysis  developments  including  new  additions, 
validation  studies,  and  conclusions.  Section  5 summarizes  the  current 
status  of  the  LASS  technique  and  recommends  future  material  modeling 
efforts  building  on  the  present  work. 

Appendix  A describes  the  LASS  analysis.  Appendix  B summarizes 
and  references  previous  results  in  the  SRI  material  properties  program. 
Appendix  C describes  the  methods  for  digital  smoothing  of  field  data, 
and  Appendixes  l)  through  F present  error  propagation  analyses. 
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2.  FIELD  TESTS 


Two  spherical  field  tests,  LASS  1 and  LASS  2,  were  performed  in 
the  Queen  15  area  of  White  Sands  Missile  Range  (WSMR)  near  the  sites 
of  the  Pre-Dice  Throw  II  events.  The  purpose  of  the  experiments  was 
to  provide  data  for  proof  testing  the  LASS  technique  at  a site  of 
interest  to  DNA.  This  section  describes  the  experiments  and  the  resulting 
data . 

We  encountered  the  usual  problems  of  field  work,  including  drill 
rig  breakdowns,  generator  failures,  instrumentation  malfunctions,  site 
lockouts,  and  stormy  weather.  However,  thanks  to  the  excellent  support 
provided  by  WES,  CERF,  WSMR,  and  DNAFC , all  phases  of  the  field  work 
were  completed  on  schedule  and  in  accordance  with  our  experimental 
requirements. 

The  experiments  were  successful.  All  the  ytterbium  stress  gages 
and  mutual  inductance  particle  velocity  gages  worked  well;  laboratory- 
quality  records  were  obtained  from  all  data  channels. 

GENERAL  REQUIREMENTS  FOR  LASS  FIELD  TESTS 

The  objective  of  a LASS  experiment  is  to  determine  the  stress-strain 
response  of  in  situ  material  under  conditions  of  dynamic  deformation. 
Significant  features  of  the  LASS  technique  are  that  the  in  situ  measure- 
ments are  made  in  relatively  undisturbed  material  and  that  the  loading 
rates  of  the  experiment  are  at  least  equal  to  the  loading  rates  of  interest 
in  stress  wave  hazard  calculations.  The  LASS  technique  is  a general  one, 
adaptable  to  cylindrical  flow,  spherical  flow,  or  plane  flow.  Both  the 
experimental  data  requirements  and  the  details  of  the  analysis  are  specific 
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to  the  flow  geometry  chosen;  it  is  crucial  to  the  validity  of  the  analysis 
that  the  actual  experimental  geometry  match  the  geometry  assumed  for  the 
analysis . 

An  ideal  spherical  LASS  experiment  would  comprise  a spherical 
explosive  charge  embedded  in  an  isotropic  medium  and  surrounded  by 
radial  stress  and  particle  velocity  gages  located  at  various  distances 
from  the  charge  center.  The  gages  are  Lagrangian;  they  move  with  the 
flow  along  particle  paths,  and  ideally,  they  do  not  perturb  the  flow. 

Since  knowledge  of  both  independent  flow  variables  (radial  stress  and 
radial  particle  velocity)  is  necessary  for  a LASS  analysis,  the  range 
of  Lagrange  coordinates  (initial  distance  from  the  charge  center)  spanned 
should  be  the  same  for  both  gage  types.  The  measured  stress  and  particle 
velocity  histories  permit  unambiguous  calculation  of  material  properties 
over  the  range  of  deformations  and  deformation  rates  spanned  by  the  gages. 
It  is  important  that  the  gages  survive  and  record  accurately  during 
stress  wave  decay,  since  the  analysis  is  valid  for  determining  unloading 
as  well  as  loading  response  of  the  material. 


FIELD  TEST  DESCRIPTION 
Site  Description 

The  Queen  15  site  of  WSMR  was  of  interest  to  DNA  because  of 

previous  work  there  including  the  two  large  Pre-Dice  Throw  II  events 

and  a CIST  experiment.  Furthermore,  extensive  field  and  laboratory 

studies  of  the  properties  of  the  site  had  been  performed  by  WES  in 

15 

support  of  the  Pre-Dice  Throw  II  events.  These  studies  had  indicated 
the  presence  of  a uniform,  fully  saturated,  clay  layer  bounded  by  the 
water  table  at  a depth  of  about  2 m and  a sandy  layer  below  at  about 
5 m.  This  layer  appeared  ideally  suited  to  a properly  designed 
spherical  I.ASS  experiment. 
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Figure  1,  a map  of  the  test  area,  shows  the  relationship  of  the 
two  LASS  experiments  to  the  CIST  and  Pre-Dice  Throw  II  events.  The 
LASS  sites  may  be  considered  as  previously  undisturbed;  peak  pressures 
from  either  of  the  Pre-Dice  Throw  II  events  or  from  the  CIST  experiment 
had  attenuated  to  less  than  a few  bars  at  the  LASS  ground  zeros  (GZs) . 

The  I.ASS  GZs  were  selected  by  DNA  Field  Command  with  SHI  concurrence. 

A basic  assumption  made  in  the  design  of  the  experiments  was  that 
the  soil  profiles  at  the  LASS  GZs  would  closely  correspond  to  the 

profiles  observed  at  test  borings  U-3,  U-7,  and  U-2  (see  Figure  1). 

To  check  on  the  validity  of  this  assumption,  we  were  alert  for  any 
anomalies  that  might  be  observable  during  (trilling  of  the  charge  holes 
and  gage  holes  for  the  LASS  experiments.  Over  the  depth  range  of 
concern,  2.4  m to  4.3  m,  the  clay  layer  appeared  to  be  uniform,  fully 
saturated,  and  of  very  low  strength.  .Results  of  the  WES  subsurface 
exploration  studies  are  given  in  reference  15. 

Experimental  Configuration 

The  general  experimental  layout  of  both  LASS  experiments  is 
illustrated  in  the  elevation  view  of  Figure  2.  The  explosive  source,  a 
center-initiated  112-kg  (256-pound)  TNT  sphere,  was  emplaced  at  a depth 
(ground  level  to  center  of  sphere)  of  3.66  m (12  feet).  Ytterbium 
stress  gages  were  grouted  into  vertical  drill  holes  at  various  distances 
from  the  charge  center.  Mutual  inductance  particle  velocity  gages  were 
grouted  into  30-degrcc  (with  respect  to  the  horizontal)  drill  holes 
radial  to  the  charge  center. 

For  practical  reasons,  various  compromises  had  to  be  made  with 
the  'ideal  experiment"  described  previously.  The  isotropy  of  the  medium 
is  perturbed  by  grout  in  the  gage  holes,  by  the  presence  of  the  cylindrical 
grouted-in  aluminum  charge  hole  liner,  and  by  the  water-saturated  sand 
surrounding  the  charge. 
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LASS  GZ  LOCATIONS,  WSMR  COORDINATES 

LASS  1 GZ:  X 490,549.78 

Y 541,370.67 

E 4,199.87 

LASS  2 GZ:  X 490,627.78 

Y 541,591.89 

E 4,201.24 

M A 46G:i  9 


FIGURE  1 MAP  OF  TEST  AREA 
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FIGURE  2 ELEVATION  VIEW  OF  EXPERIMENTAL  LAYOUT,  LASS  1 and  2 (not  to  scale) 


The  effect  of  grout  in  the  gage  holes  is  minimized  by  choosing  a grout 
that  closely  matches  the  medium,  by  spacing  the  gage  holes  to  avoid 
shadowing  of  any  gage  by  a nearby  grout  column,  and  by  restricting  the 
range  of  measurement  to  stresses  well  in  excess  of  the  grout  strength. 

The  grout  used,  DTE-2  (E) , was  recommended  by  Ralph  Bendenelli  of  WES 
for  use  in  the  saturated  clay  at  the  Pre-Dice  Throw  II  site. 

The  grouted-in  corrugated  aluminum  charge  hole  liners  and  the 
saturated  sand  surrounding  each  charge  were  required  by  the  logistics 
of  the  LASS  experiments.  To  minimize  experimental  costs,  it  was  necessary 
to  drill  the  charge  holes  and  the  gage  holes  during  the  same  period  of 
time,  about  a month  before  the  scheduled  execution  of  the  LASS  experiments. 
However,  holes  drilled  in  the  water-saturated  clay  layer  began  to  slump 
almost  as  soon  as  the  drilling  tool  was  removed.  Since  emplacement  of 
the  explosive  charges  immediately  after  drilling  would  have  been  unsafe, 
liners  were  necessary  to  keep  the  charge  holes  open  so  that  emplacement 
of  each  charge  could  be  delayed  until  the  day  before  the  scheduled  test. 
Safe  considerations  also  dictated  use  of  water-saturated  sand  to  stem 
the  charges.  Ordinary  grouts  would  not  set  up  sufficiently  in  24  hours 
to  stem  the  charge,  and  fast-setting  grouts  become  dangerously  hot. 

It  must  be  noted  that  spherical  flow  (for  a spherical  LASS  experiment) 
is  strictly  necessary  only  within  the  region  spanned  by  the  stress  and 
particle  velocity  gages.  Material  strength  effects  due  to  the  presence 
of  the  grouted  liner  are  negligible  in  the  high  stress  close-in  region, 
and  water  saturated  sand  is  a fairly  good  shock  impedance  match  to  water- 
saturated  clay.  By  centering  the  stress  gages  in  the  horizontal  plane 
passing  through  the  charge  center  and  by  placing  the  particle  velocity 
gages  at  the  relatively  shallow  angle  of  30  degrees,  we  believe  we 
avoided  the  region  in  which  the  sphericity  of  flow  was  substantially 
perturbed  by  the  presence  of  the  liner  and  by  the  small  impedance  mismatch 
between  the  sand  and  the  clay. 
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Figure  3 is  a schematic  plan  view  of  the  gage  layout  for  LASS  1. 

Note  that  the  particle  velocity  gages  actually  lie  on  a conical  surface 
30  degrees  above  the  horizontal  plane  of  the  stress  gage  centers 
(see  Figure  2).  However,  Figure  3 shows  the  true  relative  radial 
positions  of  all  gages;  the  particle  velocity  gage  positions  have  been 
rotated  into  the  plane  of  the  stress  gages.  The  gages  span  the  nominal 
range  (radial  distance  from  charge  center)  from  0.61  m (2.0  feet)  to 
1.62  m (5.3  feet).  Three  stress  gages,  120  degrees  apart,  were  placed 
at  the  same  nominal  radius  of  1.07  m (3.5  feet).  Three  particle  velocity 
gages,  all  at  the  nominal  radius  of  1.07  m,  were  placed  in  a similar 
120-degree  pattern  (rotated  60  degrees  with  respect  to  the  stress  gage 
pattern) . 

This  redundant  gage  placement  scheme  was  used  because  we  were 
concerned  about  the  accuracy  with  which  Lagrangian  gage  positions  (initial 
distance  from  the  charge  center)  could  be  determined  under  field  conditions. 
Our  concern  proved  to  be  justified;  the  saturated  clay  was  so  unstable 
that  gage  placement  and  grouting  had  to  be  done  immediately  after 
removal  of  the  tool  from  the  hole.  There  was  no  time  to  resurvey. 
Furthermore,  even  a visual  check  of  gage  position  and  hole  straightness 
was  not  possible  because  the  holes  were  over  half  full  of  water. 

We  also  encountered  problems  in  accurately  emplacing  the  TNT 
sphere  for  LASS  1.  The  charge  floated  out  of  position  as  we  were  stemming 
the  hole;  we  had  neglected  to  consider  that  the  mixture  of  sand  and 
water  could  behave  as  a dense  liquid.  We  used  a long  stick  to  push  the 
TNT  down  to  the  correct  depth,  but  we  could  no  longer  see  enough  of  the 
charge  to  tell  whether  it  was  centered  in  the  charge  hole.  Rather  than 
exhume  the  charge,  bail  out  the  hole,  and  start  over,  we  decided  to 
continue  stemming  the  hole.  We  did  not  want  to  have  to  delay  the  firing 
schedule,  and  we  were  certain  that  it  would  be  simple  to  determine  the 
true  charge  center  from  the  redundant  1.07  m particle  velocity  gage  records. 
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Figure  4 is  a schematic  plan  view  of  the  gage  layout  for  LASS  2. 


The  stress  gage  layout  is  identical  to  that  used  for  LASS  1,  but  we 
were  forced  by  circumstances  to  use  a different  particle  velocity  gage 
layout.  The  drill  rig  shifted  during  drilling  of  one  of  the  first 
nominal  30-degree  holes  and  the  resultant  25-degree  hole  was  too  far 
off-center  to  be  used  for  a close-in  gage.  We  therefore  modified  the 
layout,  dispensing  with  the  1.07-m  redundant  particle  velocity  gages. 

It  was  consequently  more  important  to  position  the  charge  accurately, 
and  we  believe  we  did  so.  We  then  blocked  the  charge  so  that  it  could 
not  float  out  of  position  while  the  hole  was  being  stemmed. 

Determination  of  Lagrange  Gage  Positions 

Accurate  knowledge  of  Lagrange  gage  positions  (initial  distance 
from  charge  center)  is  an  important  factor  in  the  LASS  analysis.  Ideally, 
the  gage  and  charge  positions  would  be  determined  by  post-installation 
surveys,  but,  as  previously  mentioned,  this  procedure  was  precluded  by 
field  conditions.  Therefore,  knowledge  of  gage  positions  depends  on 
the  accuracies  of  the  initial  survey,  of  the  drilling  of  gage  holes,  and 
of  gage  and  charge  emplacement. 

The  initial  survey  of  gage  hole  locations  was  performed  after 
drilling  and  lining  of  the  charge  holes.  The  center  of  each  liner 
was  determined  at  the  desired  charge  emplacement  depth  of  3.66  m.  This 
point  established  the  desired  GZ  (the  point  at  surface  level  directly 
above  the  charge  center)  from  which  all  gage  hole  locations  were  measured. 

For  the  vertical  stress  gage  holes,  the  drill  rig  was  spotted  with 
an  accuracy  of  about  30  mm  and  initially  aligned  to  within  \ degree  of 
vertical.  However,  because  of  drill  rig  looseness,  the  drilled  holes 
could  have  deviated  from  vertical  by  as  much  as  1 degree.  The  stress 
gages  were  attached  to  a PVC  pipe  for  installation.  A crossbar  taped 
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PLAN  VIEW  OF  GAGE  LAYOUT  FOR  LASS  2. 
Particle  velocity  gages  are  shown  as  rectangles, 
stress  gages  as  circles.  Gage  locations  to  scale. 
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to  the  pipe  provided  positive  control  of  installation  depth.  However 
the  position  of  the  gage  within  the  hole  could  not  be  controlled,  and 
the  actual  diameter  of  the  hole  at  gage  depth  was  not  known.  Our  best 
estimate  is  that  the  uncertainty  in  stress  gage  position  could  have 
been  as  large  as  ± 60  mm. 

The  30-degree  angle  of  the  particle  velocity  gage  holes  was  initially 
set  with  an  accuracy  of  1/2  degree,  the  radial  alignment  of  the  drill 
in  the  30-degree  plane  was  accurate  to  within  \ degree,  and  initial 
spotting  of  the  drill  was  accurate  to  within  about  40  mm.  PVC  pipe 
handles  marked  at  the  correct  installation  length  were  used  with  the 
particle  velocity  gages.  By  measuring  the  distance  from  the  GZ  to  the 
mark  on  the  handle,  it  was  possible  to  compensate  for  the  drill  spotting 
error.  Except  for  one  instance  in  which  a shift  of  the  drill  rig 
resulted  in  a 5-degree  deviation,  measurements  of  the  drill  stem  angle 
at  full  drill  depth  indicated  that  the  30-degree  angle  was  maintained 
with  an  accuracy  of  about  1 degree.  Allowing  for  possible  bending  of 
the  poorly  supported  portion  of  the  drill  stem  below  the  water  table, 
we  estimate  that  the  uncertainty  in  particle  velocity  gage  position 
was  no  more  than  about  20  mm. 

Because  of  their  greater  accuracy,  only  the  particle  velocity  gage 
positions  were  used  in  constructing  a time-of-arrival  versus  range 
relation.  This  relation,  together  with  stress  gage  time-of-arrival 
data,  was  used  to  obtain  more  accurate  estimates  of  stress  gage  ranges. 

For  LASS  1,  the  spread  in  times  of  arrival  for  the  three  particle 
velocity  gages  at  a nominal  range  of  1.07  m indicated  that  the  true 
charge  center  was  displaced  from  the  design  center  by  about  50  mm. 

Using  time-of-arrival  data  from  all  seven  particle  velocity  gages,  we 
iteratively  adjusted  the  charge  center  until  the  recalculated  gage  positions 
fit  a smooth  time-of-arrival  versus  range  relation. 


33 


Time-of-arrival  data,  nominal  gage  positions,  and  corrected  gage 
positions  used  in  subsequent  analyses  for  LASS  1 and  LASS  2 are  listed 


I 

> 


; 
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in  Table  1. 

Instrumentation 

The  basic  transducers  used  in  the  LASS  experiments  were  the 
ytterbium  piezoresistive  stress  gage  and  the  mutual  inductance  particle 

velocity  gage.  Both  gage  types  have  been  extensively  developed  and 

. 1-3,  5-9,  11-14 

tested  in  numerous  field  experiments.  Indeed,  the  design 

parameters  of  the  LASS  experiments  were  specifically  chosen  to  exploit 

the  stress  range  '(0.05  to  2 GPa)  over  which  both  gage  types  are  accurate 

and  have  previously  exhibited  good  survival  through  stress  release.  The 

stress  gages  were  of  the  same  designs  used  in  the  Pre-Dice  Throw  II 
14 

events,  Figure  5(b) , (c) . The  particle  velocity  gages  were  1.02  m 
long,  wound  with  PVC  insulated  solid  AWG  20  copper  wire,  Figure  5(a). 
Primary  and  secondary  windings  each  consisted  of  ten  turns.  Series 
inductances  (80  /A0  were  used  in  the  primary  circuit;  a parallel 
noninductive  50-ohm  resistor  was  also  used  in  the  primary  circuit  to 
attenuate  the  noise  pulse  that  occurs  on  gage  failure. 

The  ytterbium  stress  gages  are  powered  by  the  same  power  supplies 

2 

used  in  previous  field  experiments.  The  gage  and  its  cable  form  one 
leg  of  a bridge  circuit  powered  by  a triggered  capacitor  discharge. 

The  time  constant  of  the  discharge  exceeds  100  ms;  voltage  across  the 
bridge  is  essentially  constant  over  the  few  milliseconds  spanning  gage 
recording  times.  The  gage  changes  resistance  upon  arrival  of  the  stress 
wave,  and  the  resultant  bridge  imbalance  is  recorded  by  an  oscilloscope 
camera,  backed  up  by  a tape  recorder.  Before  the  experiment,  the  system 
response  is  calibrated  by  series  insertion  of  precision  resistors  spanning 
the  range  of  expected  gage  resistance  changes.  These  calibrations  are 
used  to  convert  the  observed  voltage-time  records  into  relative 
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Table  1 


RANGE  CORRECTIONS  FOR  LASS  1 AND  LASS  2 


Nominal 

Time  of 

Corrected 

Gage 

Range 

Arrival 

Range 

No. 

(m) 

Cite) 

(m) 

/ 

' 1 

0.610 

144 

0.526 

1 

2 

0.671 

204 

0.692 

LASS  1 1 

3 

0.732 

189 

0.655 

Stress  gages  \ 

4 

1.07 

334 

0.987 

1 

5 

1.07 

313 

0.942 

[ 

6 

1.07 

343 

0.999 

7 

1.615 

612 

1.580 

( 12 

0.610 

175 

0.610 

| 

13 

0.732 

242 

0.701 

LASS  1 ! 

10 

0 . 884 

290 

0.930 

Particle  velocity  ^ 

, 14 

1.07 

356 

1.02 

11 

1.07 

388 

1.10 

gages  | 

8 

1.07 

404 

1.10 

\ 9 

1.615 

624 

1.65 

0.610 

150 

0.540 

I 2 

0.671 

168 

0.611 

LASS  2 

) 3 

0.732 

176 

0.631 

Stress  gages 

1 

1.07 

308 

0.950 

1 5 

1.07 

309 

0.959 

6 

1.07 

318 

0.976 

V 7 

1.615 

572 

1.510 

i 

( 12 

0.610 

169 

0.610 

| 

13 

0.732 

219 

0.732 

LASS  2 

10 

0 . 884 

267 

0.884 

Particle  velocity  i 

i 14 

1.07 

366 

1.07 

gages  i 

11 

1.31 

474 

1.31 

1 

( « 

1.49 

596 

1.49 

' 9 

1.615 

616 

1.615 
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resistance-time  data,  with  appropriate  corrections  for  cable  resistance. 

Stress-time  data,  for  both  loading  and  unloading,  are  obtained  from  the 

9 

ytterbium  calibrations  of  Ginsberg  et  al. 

The  mutual  inductance  gages  are  also  powered  by  capacitor 

2 

discharge  power  supplies,  as  described  in  a previous  report.  In  the 
present  experiments,  the  discharge  time  constant  was  only  about  5 ms, 
i.e.,  only  an  order  of  magnitude  above  gage  recording  time.  Consequently, 
gage  currents  decayed  by  about  10%  over  the  times  spanned  by  the  particle 
velocity  gage  records. 

Preshot  calibrations  of  each  gage  system  indicated  that  the  current- 
versus-time  relationship  was  reproducible  to  within  about  1%.  The 
experimental  data  are  voltage-time  histories,  recorded  by  oscilloscope 
cameras  backed  up  by  tape  recorders.  These  voltage  histories  are  conver- 
ted to  particle  velocity  histories  by  the  following  equation: 


T v M . 

v<t>  = l(t>  -«<t> 

T C 

where  R^  is  the  terminating  resistor  on  the  recording  end  of  the  signal 
cable,  R^  is  the  resistance  of  the  signal  cable,  M is  the  mutual 
inductance  of  the  gage,  and  X is  the  length  of  the  gage  coils.  V(t), 

I(t),  and  u(t)  are  the  voltage,  current,  and  particle  velocity, 
respectively.  V(t)  is  measured  with  respect  to  the  zero-particle-velocity 
baseline  measured  in  preshot  calibrations.  Because  of  the  time-varying 
current  and  the  inductive  nature  of  the  gage,  the  baseline  itself  is  a 
time-varying  voltage. 
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Instrumentation  Improvements  for  LASS  Experiments 


One  of  the  major  objectives  of  the  field  experiments  was  to  obtain 
records  of  sufficiently  good  quality  to  permit  taking  full  advantage  of 
the  LASS  te- Unique.  Problems  encountered  in  previous  field  experiments 
included  premature  gage  failure  (i.e.,  before  complete  unloading)  at 
stresses  above  about  1 GPa , severe  electrical  noise  pickup,  and  poor 
signal-to-noise  ratio  for  mutual  inductance  particle  velocity  gages  at 
low  particle  velocities.  Electrical  crosstalk  between  gages  was  also 
a problem.  The  improvements  described  below  allowed  us  to  obtain 
laboratory-quality  gage  records  from  both  LASS  experiments. 

Premature  gage  failure  has  been  traced  to  a variety  of  causes 

including  shear  failure  at  the  gage  package-cable  junction,  voids  in 

the  grout  column,  and  cracking  of  the  grout  column  (due  to  excessive 

moisture  loss) . The  shear  failure  problem  appears  to  have  been  solved 

by  simple  modifications  in  the  construction  of  the  gage  package,  shown 

in  Figure  5(b) . The  RG-8  lead  exits  the  package  at  a 30-degree  angle 

(rather  than  in  the  plane  of  the  package)  and  loops  around  to  join  the 

PVC  pipe  used  as  an  installation  handle.  This  package  design  was  also 

used  with  good  results  in  the  100-ton  TNT  and  120-ton  ANFO  Pre-Dice 
. 14 

Throw  II  events.  Voids  in  the  grout  column  were  avoided  by 
careful  grouting  from  the  bottom  of  the  hole,  and  the  choice  of  a wet 
site  circumvented  the  grout  cracking  problem.  Figure  6 is  a reproduction 
of  a LASS  1 stress  gage  record  showing  recording  duration  extending 
through  complete  unloading  from  a peak  stress  of  1.71  GPa. 

The  major  instrumentation  improvement  in  the  WSMR  LASS  experiments 
was  the 'modification  of  the  mutual  inductance  particle  velocity  gage 
power  supplies  for  remote  operation.  This  allowed  us  to  place  the 
power  supplies  near  the  experiment  rather  than  in  the  instrumentation 
van,  as  had  been  previous  practice.  The  power  supplies  were 
placed  in  a plywood  shelter,  protected  by  sandbags,  about  50  m from  the 
LASS  2 GZ  (Figure  1) . 
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Close-in  power  supply  placement  effectively  eliminates  the  problem 
of  noise  pickup  caused  by  inductive  coupling  between  power  supply  leads 
and  signal  leads.  The  IR  drop  in  the  power  leads  is  reduced,  and  since 
the  gage  signal  depends  linearly  on  gage  current,  the  signal-to-noise 
ratio  at  low  particle  velocities  is  substantially  improved.  An  additional 
advantage  is  a substantial  reduction  in  cable  costs.  In  previous  field 
tests  it  was  necessary  to  use  high-cost  low-resistance  cables  (RG-331) 
for  power  supply  leads  whenever  the  run  was  over  about  300  m.  Even  so, 
typical  gage  currents  were  limited  to  50-70  amperes  by  cable  resistance. 

In  contrast,  gage  currents  in  the  present  experiments  were  about  175 
amperes.  Figure  .7  is  a reproduction  of  a LASS  1 MIPV  gage  record,  showing 
good  signal-to-noise  ratio  even  at  a particle  velocity  of  about  15  m/s. 

Major  sources  of  electrical  crosstalk  between  gages  are  ground 
loops  and  inductive  coupling  between  cables.  The  ground  loop  problem 
was  minimized  by  floating  all  the  power  supplies  and  the  instrumentation 
van.  A ground  loop  can  be  created  when  a gage  shorts  to  ground  as  a 
result  of  shock  wave  damage.  If  this  short  occurs  early  with  one  of 
the  close-in  gages,  the  records  of  the  gages  further  out  can  be  affected. 
With  a floating  instrumentation  van,  two  close-in  gages  must  short  to 
ground  before  a ground  loop  can  form. 

Inductive  coupling  between  cables  was  negligible  in  LASS  1.  We  had 
time  to  achieve  a very  clean  cable  layout.  Each  signal  cable  was 
separated  from  parallel  adjacent  cables  by  at  least  15  cm  over  the 
entire  run  from  shot  to  instrumentation  van,  and  all  cable  crossings 
were  at  right  angles.  The  IASS  1 cables  were  reused  for  LASS  2,  and  we 
had  to  forego  a clean  cable  layout  in  order  to  meet  the  firing  schedule. 
Although  the  LASS  2 records  are  much  better  than  previous  field  records, 
they  are  inferior  to  the  LASS  1 records. 
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FIELD  TEST  DATA  AND  DISCUSSION 


Composite  experimental  data  for  LASS  1 and  LASS  2 are  presented 
in  Figures  8 through  11,  and  the  individual  gage  profiles  are  presented 
in  Figures  12  through  15.  The  original  gage  records  were  digitized  on 
a Telerecordex.  Simple  computer  programs  apply  the  calibration  data  to 
the  digitized  records  to  convert  the  measured  voltage  histories  to  stress 
and  particle  velocity  histories. 

To  prepare  the  data  for  further  reduction  and  to  prepare  Figures 
8-11,  we  redigitized  the  gage  profiles  at  uniform  time  steps  of  2.5  Jls 
and  removed  obvious  noise  glitches.  The  solid  lines  of  Figures  12 
through  15  show  the  original  profiles,  and  the  dashed  lines  show  the 
profiles  after  noise  removal.  Noise  glitches  can  be  unambiguously 
identified  by  their  simultaneous  occurrence  on  two  or  more  records.  In 
general,  a noise  burst  at  late  time  coincides  with  the  failure  of  one 
of  the  close-in  gages. 

The  results  of  the  WSMR  LASS  experiments  demonstrate  that  it  is 
possible,  under  field  conditions,  to  obtain  good  quality,  long  duration, 
stress  and  particle  velocity  records  over  the  stress  range  between  0.05 
and  2 GPa.  However,  all  the  problems  of  field  measurements  have  not 
been  solved.  Further  attention  should  be  given  to  the  problem  of  gage 
emplacement  in  dry  soils,  in  which  drying  out  and  cracking  of  grout 
columns  can  lead  to  early  gage  failure.  Stress  gage  systems  suitable 
for  use  in  the  range  below  about  0.05  GPa  have  not  vet  been  adequately 
developed,  although  standard  accelerometers  can  be  used  for  particle  velocity 
measurements  in  this  range. 

We  were  able  to  minimize  the  electrical  noise  problem  because  we 
had  good  control  over  the  instrumentation  and  cabling.  However,  in  the 
typical  large-scale  field  experiment,  instrumentation  responsibilities 
are  divided  among  several  agencies,  and  mutual  interference  is  the  rule. 
Nevertheless,  current  field  recording  systems  for  both  stress  gages  and 
particle  velocity  gages  could  be  improved  to  reduce  their  sensitivity  to 
resulting  extraneous  electrical  fields. 


■12 


TIME,  n s 

MA-4663-2 

FIGURE  9 PARTICLE  VELOCITY  HISTORIES,  LASS  1. 

Gages  located  0.610,  0.701 , 0.930,  1 .02,  1.10,  and  1 .65  m from  charge  center. 
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FIGURE  10  RADIAL  STRESS  HISTORIES,  LASS  2. 

Gages  located  0.540,  0.61 1 , 0.631 , 0.950,  0.959,  0.976,  and  1 .51  m from  charge  center. 
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FIGURE  13  INDIVIDUAL  PARTICLE  VELOCITY  GAGE  RECORDS,  LASS  1 
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FIGURE  15  INDIVIDUAL  PARTICLE  VELOCITY  GAGE  RECORDS,  LASS  2 (Concluded) 


3.  LASS  ANALYSIS  OF  FIELD  TESTS 


The  piecewise  linear  form  of  LASS  (Appendix  A)  was  used  to  analyze 
data  from  the  two  spherical  field  tests,  LASS  1 and  LASS  2.  The  objectives 
were  to  proof  test  LASS  using  field  data  from  a site  of  interest  to  DNA 
and  to  obtain  dynamic  in  situ  material  property  information  for  the 
test  material,  wet  clay  from  the  Pre-Dice  Throw  II  site.  Both  objectives 
were  satisfied. 

SMOOTHING 

The  first  step  in  this  form  of  LASS  is  to  presmooth  the  records  to 
reduce  random  high  frequency  noise.  The  field  test  records  were  smoothed 
with  the  new  modified  low-pass  filter  technique  described  in  Section  4. 
Specific  smoothing  parameters  are  given  in  Appendix  C.  The  resulting 
smoothed  records  are  shown  in  Figures  16  through  19.  These  may  be  com- 
pared with  the  original  data,  Figures  8 through  11,  to  judge  the  rela- 
tively small  amount  of  smoothing  imposed  (significant  features  are  not 
lost) . 

LASS  1 ANALYSIS  AND  RESULTS 

The  smoothed  stress  profiles  from  LASS  1,  Figure  16,  were  used  to 
* 

generate  a five-sheet  gage  flow  surface  for  the  piecewise  linear  analysis. 

Five  characteristic  features,  or  knots,  in  addition  to  the  foot  of  the 

wave  were  identified  on  each  record.  Sequentially  from  the  foot  these 

were:  the  peak  CT  , feature  A,  feature  B,  feature  C, 

max 

♦ 

This  terminology  refers  to  the  number  of  piecewise  linear  sections,  or 
sheets  used  to  characterize  the  gage  flow  surface.  Figure  A4  in  Appendix 
A shows  a four-sheet  gage  flow  surface  with  four  piecewise  linear  gage 
traces  on  it  crossing  each  sheet. 
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FIGURE  16  SMOOTHED  RADIAL  STRESS  HISTORIES,  LASS  1. 

Gages  located  0.526,  0.655,  0.692,  0.942,  0.987,  0.999,  and  1.58  m from  charge  center. 
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FIGURE  17  SMOOTHED  PARTICLE  VELOCITY  HISTORIES,  LASS  1. 

Gages  located  0,610,  0.701 , 0.930,  1 .02,  1.10,  and  1 .65  m from  charge  center. 
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FIGURE  19  SMOOTHED  PARTICLE  VELOCITY  HISTORIES,  LASS  2. 

Gages  located  0.610,0.732,  0.884,  1.07.  1.31,  1.49,  and  1.62  m 
from  charge  center 


and  0.1  CT  . These  features  are  illustrated  on  an  individual  record 
max 

in  Figure  20.  These  particular  five  knots  were  selected  because  they 
adequately  represent  the  stress  gage  flow  surface  and  because  they 
result  in  a satisfactory  representation  of  the  particle  velocity  gage 
flow  surface,  as  discussed  later. 

Next,  the  relative  time  (time  measured  from  the  foot  of  the  wave 

at  that  Lagrange  coordinate)  and  the  amplitude  for  each  knot  were 

smoothly  connected  in  the  spatial  direction  by  least-squares  fitting  to 

functions  of  Lagrange  coordinates,  h.  Specifically,  the  time  of  arrival 

* 

of  the  foot  of  the  wave  was  fit  to  a linear  function  of  Lagrange  coor- 

t 

dinate  h (initial  particle  position)  and  the  rest  of  the  features  were 
fit  to  linear  or  quadratic  exponential  functions  of  h.  The  specific 
analytic  functions  and  parameters  used  for  the  stress  link  lines  are 
given  in  Table  2.  The  resulting  link  lines  and  data  points  are  shown 
in  Figures  21  and  22. 

The  smoothed  LASS  1 particle  velocity  profiles,  Figure  17,  were 
also  characterized  by  five  knots  in  addition  to  the  foot  of  the  wave  to 


Since  the  time  of  first  motion  in  stress  and  particle  velocity  is  the 
same,  both  stress  and  particle  velocity  data  are  generally  used  to  fit 
this  feature.  In  the  present  case,  however,  because  of  the  soupy  nature 
of  the  site  material,  the  locations,  h,  of  the  particle  velocity  gages 
were  known  more  accurately  than  those  of  the  stress  gages,  as  discussed 
in  Section  2.  Consequently,  only  particle  velocity  gage  data  were  used 
to  fit  the  foot  link  line. 

t 

Since  the  relative  times  at  which  the  stress  and  particle  velocity 
maxima  occurred  appeared  to  scatter  about  the  same  curve  withir  experi- 
mental error,  data  from  both  gage  types  were  used  to  fit  this  feature 
and  the  resulting  link  line  was  used  to  represent  the  peak  in  both 
gage  flow  surfaces.  When  the  data  scatter  about  a single  curve,  fitting 
them  together  improves  accuracy  as  discussed  in  reference  1. 
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LASS  1 STRESS  GAGE  FLOW  SURFACE  LINK  LINES 
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MA  4663  25 

FIGURE  21  TIME  LINK  LINES  FOR  LASS  1 STRESS  GAGE  FLOW  SURFACE. 
Equations  for  these  link  lines  are  given  In  Table  2. 
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generate  a five-sheet  gage  flow  surface.  However,  these  knots  were  not 
selected  independently  of  the  stress  gage  features.  Instead,  at  each 
particle  velocity  gage  location,  h,  the  time  of  each  of  the  five  fitted 
stress  gage  flow  knots  was  determined  from  Figure  21(b).  At  each  of 
these  five  times  on  each  smoothed  particle  velocity  record,  the  amplitude 
was  determined.  These  are  the  five  knots  used  to  represent  the  particle 
velocity  records;  they  are  illustrated  on  one  of  the  individual  particle 
velocity  records  in  Figure  23.  This  procedure,  developed  in  the  previous 
year  to  improve  the  stress  difference  calculation,  ensures  compatibility 
of  stress  and  particle  velocity  derivatives  at  joints  in  the  respective 
gage  flow  surfaces.  With  this  procedure,  the  relative  time  link  lines 
are  the  same  ones  used  for  the  stress  gages  (Table  2 and  Figure  21). 

The  link  line  functions  used  to  fit  the  particle  velocity  amplitudes 
are  given  in  Table  3 and  shown  in  Figure  24.  Stress  and  particle  velocity 
records  generated  from  the  link  lines  in  Tables  2 and  3 are  given  in 
Figures  25  and  26  for  comparison  to  the  input  profiles  in  Figures  16  and 
17. 

The  link  lines  in  Tables  2 and  3 and  an  initial  density  value  of 
3* 

1.85  g/cm  were  used  in  the  LASS  subroutine  STPV  to  calculate  the  stress  and 
strain  histories  of  different  particles  during  the  LASS  1 event.  From 
these,  any  desired  stress-strain,  stress-stress,  or  strain-strain  cross- 
plot can  be  constructed.  The  results  are  given  in  Figures  27  through  29. 

Figure  27  shows  radial  stress  versus  volume  strain  at  five  particle 
locations  within  the  instrumented  region  of  the  flow.  For  weak  site 
material  at  high  stresses  (the  present  case),  radial  stress  is  expected 
to  be  very  close  to  mean  stress,  so  Figure  27  is  equivalent  to  a mean 
stress-volume  strain  plot.  This  is  often  the  most  important  stress- 
strain  information  for  modeling  purposes.  Figure  28  gives  three  other 
♦ 

Private  communication,  A.  F,.  Jackson,  WES. 
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FIGURE  23  FEATURES  USED  TO  CHARACTERIZE  LASS  1 PARTICLE  VELOCITY  RECORDS 
FOR  PIECEWISE  LINEAR  LAGRANGE  ANALYSIS. 

These  features  are  derived  from  the  stress  records  as  discussed  in  the  text.  The  featuies 
are  shown  on  the  particle  velocity  record  from  0.610  m. 


h in  meters, 
b 

Features  are  determined  from  the  stress  gage  flow  surface,  as  discussed 
in  text. 
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FIGURE  24  AMPLITUDE  LINK  LINES  FOR  LASS  1 PARTICLE  VELOCITY  GAGE  FLOWSURFACE 
Equations  for  these  link  lines  aie  given  in  Table  3. 


STRESS,  GPa 


FIGURE  25  PIECEWISE  LINEAR  GAGE  FLOW  SURFACE  REPRESENTATIONS  OF  LASS  1 STRESS 
HISTORIES  AT  LAGRANGE  COORDINATES  0.6,  0.8,  1.05,  1.3,  and  1 6 m (measured 
from  charge  center) 
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FIGURE  27  RADIAL  STRESS  VERSUS  VOLUME  STRAIN.  LASS  1 
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FIGURE  28  RADIAL  STRAIN  VERSUS  TANGENTIAL  STRAIN,  VOLUME  STRAIN  VERSUS 
PRINCIPAL  STRAIN  DIFFERENCE,  AND  RADIAL  STRESS  VERSUS  RADIAL 
STRAIN,  LASS  1,  AT  LAGRANGE  COORDINATES  0.6,  0.8,  1.05,  1.3,  and  1.6  m 
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stress-strain  or  strain-strain  cross  plots  also  of  interest  for  modeling, 
namely,  radial  strain/tangential  strain,  volume  strain/strain  difference, 
and  radial  stress/ radial  strain.  Figure  29  gives  the  individual  stresses 
and  strains  as  function  of  time  for  the  particle  at  the  middle  of  the 
instrumented  interval.  These  trajectories  are  available  for  any  particle; 
those  for  the  middle  of  the  analyzed  interval,  1.05  m,  have  been  repro- 
duced in  Figure  29.  From  these  histories  and  the  radial  stress-time 
profiles  in  Figure  25,  any  stress-strain  cross-plots  such  as  those  in 
Figures  27  and  28  can  be  constructed.  The  results  in  Figures  27  through 
29  are  discussed  below.  Note  that  in  all  LASS  analyses,  results  reported 
at  the  middle  of  the  h-interval  should  be  considered  most  accurate  because 
numerical  approximations  are  best  there. 

LASS  2 ANALYSIS  AND  RESULTS 

The  IASS  2 analysis  approach,  except  for  details,  is  the  same  as 
that  used  for  LASS  1 and  the  present  discussion  parallels  that  for  LASS 
1.  The  smoothed  LASS  2 radial  stress  profiles,  Figure  18,  were  char- 
acterized by  six  features  in  addition  to  the  foot  of  the  wave:  the 

peak  a , A,  B,  C,  0.2  a , and  0.1  cr  . These  features  are 
max  max  max 

illustrated  on  an  individual  record  in  Figure  30.  The  amplitude  and 
relative  time  link  lines  used  to  connect  these  knots  are  given  in  Table  4 
and  shown  graphically  in  Figures  31  and  32.  Note  that  as  in  LASS  1 the 
time  of  arrival  (real  time,  not  relative  time)  of  the  foot  of  the  wave 
was  fit  to  a linear  function  using  the  particle  velocity  gage  data  only. 
Also,  as  in  LASS  1,  the  relative  time  of  occurrence  of  stress  and  particle 
velocity  peaks  were  fit  together  by  a single  function  since  they  scatter 
about  a single  curve. 

The  smoothed  LASS  2 particle  velocity  profiles,  Figure  19,  were 
characterized,  as  in  IASS  1,  by  their  amplitudes  at  the  relative  times 
of  the  fitted  knots  in  stress.  The  resulting  knots  on  an  individual 
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FIGURE  30  FEATURES  USED  TO  CHARACTERIZE  LASS  2 STRESS  RECORDS  FOR  PIECEWISE 
LINEAR  LAGRANGE  ANALYSIS. 

The  six  features  are  shown  on  the  stress  record  from  0.540  m. 
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LAGRANGE  COORDINATE  (h),  m 
(a)  Foot  of  wave. 


(b)  Six  characteristic  features 

MA  4G63  3b 

FIGURE  31  TIME  LINK  LINES  FOR  LASS  2 STRESS  GAGE  FLOW  SURFACE. 
Equations  for  these  link  lines  are  given  in  Table  4 
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particle  velocity  record  are  shown  in  Figure  33.  Because  of  the  pro- 
cedure used,  the  relative  time-h  link  lines  are  the  same  for  the  stress 
and  particle  velocity  records;  this  was  also  the  case  in  the  LASS  1 
analysis.  The  link  lines  fit  to  the  particle  velocity  amplitudes  at 
the  knots  are  given  in  Table  5 and  shown  in  Figure  34.  The  stress  and 
particle  velocity  record  representations  resulting  from  this  analysis 
are  shown  in  Figures  35  and  36. 

The  link  lines  in  Tables  4 and  5 and  the  initial  density  value  of 
3 * 

1.85  g/cm  were  input  into  the  LASS  stress-  and  strain-calculating 
subroutines  to  complete  LASS  2 analysis.  The  results  are  presented  in 
Figures  37  through  39.  Because  the  quality  of  the  data  from  LASS  2 was 
lower  than  that  from  LASS  1,  it  was  necessary  to  reduce  the  h-interval 
over  which  stresses  and  strains  were  calculated.  Consequently,  when 
comparing  the  LASS  1 and  2 results,  one  must  take  care  that  the  results 
apply  to  the  same  particle  position,  h,  in  the  flow. 

DISCUSSION 

The  results  from  LASS  2 are  qualitatively  the  same  as  those  from 
LASS  1.  The  following  observations  regarding  deformation,  material 
response,  and  accuracy  are  based  on  data  from  LASS  1,  but  they  apply 
equally  to  the  results  of  LASS  2. 

Deformation 

The  idealized  deformation  of  the  site  material  in  these  spherical 

j 

experiments  is  nearly  pure  radial  strain  in  compression  followed  by 

nearly  constant  radial  strain  during  release  (see,  for  example,  Figure  28(a)). 

Note  that  these  conclusions  can  also  be  inferred  from  other  plots  of  the 

• 

Private  communication,  A.  E.  Jackson,  WES. 
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PARTICLE  VELO 


100  200  300  400  500  600 

TIME,  p s 

MA  4663  37 

FIGURE  33  FEATURES  USED  TO  CHARACTERIZE  LASS  2 PARTICLE  VELOCITY  RECORDS 
FOR  PIECEWISE  LINEAR  LAGRANGE  ANALYSIS. 

These  features  are  derived  from  the  stress  records  as  discussed  in  the  text. 

The  features  are  shown  on  the  particle  velocity  record  fiom  0.732  m. 


Table  5 


LASS  2 PARTICLE  VELOCITY  GAGE  FLOW  SURFACE  LINK  LINE3 


Amplitude  Link  Line 
Feature*3  x p,s/'mm 

Foot 

Peak,  u exp(-0 .9186  - 1.315b) 

III  cl  A 


°-2 

°-l  W 


exp(-l .326  - 1 . 099h) 
exp(-l .346  - 1.286h) 
exp(-1.835  - 1.019h) 
exp (-0 . 9548  - 2.180h) 
cxp(-0 . 9844  - 2.42910 


Relative  Time  Link  Line 
x ps-1 


Same  as  coefficients 
for  Ci  Table  4 


h in  meters . 

Features  are  determined  from  the  stress  gage  flow  surface  as 
discussed  in  text. 
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M A -4663  39 

FIGURE  35  PIECEWISE  LINEAR  GAGE  FLOW  SURFACE  REPRESENTATIONS  OF  LASS  2 STRESS 
HISTORIES  AT  LAGRANGE  COORDINATES  0.8,  0.9,  1.05,  1.2,  AND  1.3  m 
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PARTICLE  VELOCITY,  ms 


TIME,  ns 

MA-4663-40 

FIGURE  36  PIECEWISE  LINEAR  GAGE  FLOW  SURFACE  REPRESENTATION  OF  LASS  2 
PARTICLE  VELOCITY  HISTORIES  AT  LAGRANGE  COORDINATES  0.8,  0.9, 
1.05,  1.2,  AND  1.3  m 


RADIAL  STRESS,  GPa 


MA  4663  41 


FIGURE  37 


RADIAL  STRESS  VERSUS  VOLUME  STRAIN,  LASS  2,  AT  LAGRANGE 
COORDINATES  0.8,  0.9,  1.05,  1.2,  AND  1.3m 


RADIAL  STRAIN, 


0 0.02  0.04  0.06 

RADIAL  STRAIN,  er 


(c)  Radial  stress  versus  radial  strain 


MA  4663  42 


FIGURE  38  RADIAL  STRAIN  VERSUS  TANGENTIAL  STRAIN,  VOLUME  STRAIN  VERSUS 
PRINCIPAL  STRAIN  DIFFERENCE,  AND  RADIAL  STRESS  VERSUS  RADIAL 
STRAIN,  LASS  2,  AT  LAGRANGE  COORDINATES  0.8,  0.9,  1.05,  1.2,  AND  1.3  m 


H(> 


VOLUME  STRAIN,  tv  RADIAL  STRAIN, 


0.05 
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(a)  Radial  strain  versus  time 
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(b)  Tangential  strain  versus  time 


(c)  Volume  strain  versus  time 
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(d)  Stress  difference  versus  time 


M A -4663-43 

FIGURE  39  RADIAL  STRAIN,  TANGENTIAL  STRAIN,  VOLUME  STRAIN,  AND  STRESS 
DIFFERENCE  VERSUS  RELATIVE  TIME  AT  LAGRANGE  COORDINATE 
h = 1.05  m,  LASS  2 


X7 


r i 


LASS  results,  such  as  the  individual  strain  histories  in  Figure  29  and 
radial  stress-radial  strain  cross-plot  as  in  Figure  28(c).  Although 
the  idealized  deformation  is  one-dimensional  (radial)  in  compression, 
it  is  not  uniaxial  unless  the  radius  of  curvature  of  the  wave  front  at 
the  particle  in  question  is  effectively  infinite. 

The  fact  that  the  radial  strain  is  nearly  constant  in  release  can 
also  be  determined  directly  from  the  composite  plots  of  the  raw  particle 
velocity  records,  Figure  9.  The  release  profiles  for  each  gage  lie 
approximately  along  a single  line,  meaning  that  at  any  instant  the  radial 
particle  velocities  of  all  particles  undergoing  release  are  equal. 
Therefore  there  is  no  change  in  the  radial  dimension  (thickness)  of  a 
small  volume  element  and  hence  no  change  in  radial  strain  during  this 
portion  of  the  flow. 

Material  Response 

For  the  purpose  of  discussing  material  response,  we  have  reproduced 
the  radial  stress-volume  strain  plots  of  Figure  27  and  added  an  "equili- 
brium compression  curve,"  the  dark  line  in  Figure  40.  This  is  the  locus 
of  end  or  peak  compression  states  achieved  at  each  particle  during  load- 
ing. In  the  limits  of  instantaneous  loading  and  uniaxial  strain,  this 
curve  appi'oaches  the  familiar  Hugoniot.  Remember  that  the  individual 
trajectories  calculated  by  LASS  represent  the  actual  loci  of  states 
traversed  by  individual  particles  during  the  dynamic  load-release  cycle; 
the  particles  do  not  load  up  along  the  "equilibrium  compression  curve" 
in  high  rate  loading. 

Note  first  in  Figure  40  that  the  equilibrium  compression-dynamic 
release  paths  are  nearly  reversible  until  release  to  low  stresses;  i.e., 
the  equilibrium  behavior  of  the  material  is  not  hysteretic.  The  loading 
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MA -4663  44 


FIGURE  40  RADIAL  STRESS  VERSUS  VOLUME  STRAIN,  LASS  1,  FROM  FIGURE  27  WITH 
EQUILIBRIUM  COMPRESSION  CURVE  ADDED  (DARK  LINE) 
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strain  rate  varies  by  approximately  a factor  of  10  among  the  trajectories 
4-1  3-1 

from  about  10  s at  0.6  m to  about  10  s at  1.6  m.  Consequently  to  the 
extent  that  the  reversibility  observation  is  valid, and  over  this  strain 
rate  range,  the  LASS  1 material  is  rate  insensitive.  The  data  in  Figure 
40  indicate  not  only  that  the  site  material  is  not  hysteretic  at  high 
stresses,  but  that  at  low  stresses  it  has  not  been  irreversibly  compacted. 
In  fact,  the  trajectories  indicate  some  bulking  at  zero  stress.  Finally, 
note  that  the  relatively  high-stress  trajectories  show  no  indication  of 
any  site  material  strength  or  of  initial  high  compressibility.  The 
individual  stress  profiles  also  showed  no  precursor  or  other  evidence  of 
strength. 

In  summary  the  results  in  Figure  40  suggest  (1)  reversible  equili- 
brium compression-release  behavior,  (2)  no  compaction  either  initially 
or  upon  release,  and  (3)  no  detectable  strength  at  these  stresses.  This 
behavior  is  very  reasonable  for  clay  material  below  the  water  table. 
Quantitative  parameters  can  be  evaluated  from  the  LASS  plots  once  a 
specific  model  is  selected.  Comparison  of  the  LASS  results  with  data 
determined  by  other  tests  for  the  same  material  is  given  in  the  Conclu- 
sions. 

The  radial  stress  level  of  the  LASS  field  tests  is  several  orders 

of  magnitude  greater  than  the  stress  difference,  $,  that  wet  clay  can 

be  expected  to  support.  Previous  work  in  the  LASS  program  as  well  as 

the  results  presented  in  Section  4 of  this  report  show  that  in  this 

situation  (i.e.,  when  c » $) , $ cannot  be  calculated  with  precision 

r 

by  Lagrange  analysis  because  the  flow  is  insensitive  to  § under  these 
conditions.  This  is  also  apparent  by  noting  that  $ appears  only  in  the 
radial  momentum  equation  (Appendix  A)  where  it  is  expressed  as  the 
difference  of  two  derivatives  that  must  be  evaluated  from  field  data. 

When  these  derivatives  are  large  and  § is  small  (the  present  case), 


small  errors  in  the  field  data  or  its  representation  cause  large  errors 
in  the  $ calculation. 
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In  spite  of  the  large  expected  errors,  we  did  calculate  $ from  the 
LASS  1 field  data.  The  resulting  §-t  history  for  the  particle  at  the 
center  of  the  interval  is  shown  in  Figure  29.  The  initial  large  positive 
excursion  of  about  20  y,s  duration  occurs  during  compression  when  the 
flow  derivatives  and  experimental  errors  in  determining  them  are  largest. 
During  release,  the  derivatives  and  the  errors  in  estimating  them  are 
significantly  smaller.  In  this  region  the  § estimate  shows  much  smaller 
excursions  (maximum  about  0.3  GPa,  average  about  0.15).  These  are 
estimated  to  be  comparable  to  or  less  than  the  error  propagated  in  the 
calculation  (see  following  error  discussion).  Consequently,  after  the 
first  20  y,s,  the  results  suggest  that  § is  bounded  by  about  + 0.1,  - 0.3 
GPa.  Therefore,  to  the  precision  and  on  the  scale  of  these  high  level 
tests,  the  result  § equal  zero  is  permitted.  This  is  consistent  with 
the  behavior  expected  for  a wet  clay  for  which  § should  be  on  the  order 
of  10-3  GPa  or  less. 


Errors 

Three  sources  of  error  in  the  LASS  results  have  been  identified 
and  considered.  These  are  (1)  deviations  from  the  assumed  flow  symmetry 
(spherical  in  this  case),  (2)  errors  in  taking  the  field  measurements 
and  reducing  them  to  the  form  in  Figures  8 through  11  where  they  are 
ready  for  LASS  analysis,  and  (3)  numerical  or  modeling  errors  related  to 
the  representation  of  the  records  (piecewise  linear  in  this  case)  and  to 
the  numerical  approximations  in  the  LASS  analysis.  After  examining  the 
field  records,  especially  those  from  gages  located  at  similar  ranges 
but  along  different  azimuths,  we  consider  the  first  error  source 
(deviations  from  spherical  symmetry)  to  be  insignificant  compared 
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with  the  other  sources  tor  both  LASS  field  tests.  The  magnitudes  of 
these  other  errors  are  discussed  below. 

To  estimate  the  effect  of  input  errors  on  the  calculations,  we 
used  results  of  the  global  Lagrange  analysis  and  error  propagation 
feasibility  demonstration  performed  on  the  LASS  1 data  and  described  in 
Section  4.  This  is  the  first  time  any  quantitative  error  propagation 
analysis  has  been  attempted  for  LASS  calculations,  and  the  resulting 
estimates  are  meant  to  indicate  only  analysis  sensitivity  and  order  of 
magnitude  of  the  errors.  The  approach  is  not  rigorous  since  different 
analyses  (global  and  piecewise  linear)  have  been  used  to  estimate  the 
errors  and  to  calculate  the  results,  since  the  global  analysis  itself 
was  quite  crude  (see  Section  4),  and  since  the  errors  in  the  input  are 
modeled  in  the  simplest  way. 

The  LASS  input  consists  of  the  presmoothed  stress  and  particle 
velocity  records  and  the  gage  locations.  The  dominant  sources  of  errors 
in  these  data  are  taken  to  be  (1)  amplitude  errors  in  the  smoothed  gage 
records  resulting  from  gage  package  variations,  calibration  uncertainties, 
and  so  on,  and  (2)  uncertainties  in  the  gage  ranges  from  the  source. 

For  the  error  propagation  analysis,  we  took  the  standard  deviation  in 
smoothed  velocity  records  to  be  5 m s-1,  in  smoothed  stress  records  to 
be  0.013  GPa , and  in  the  gage  locations  to  be  0.038  m.  The  resulting 
standard  deviations  in  the  calculated  strains  and  stresses  are  shown  in 
Section  4.  They  may  be  taken  as  about  15%  of  the  peak  strains  and  as 
about  0.5  GPa  for  the  stress  difference  or  about  40  times  the  assumed 
error  in  the  radial  stress  measurements.  This  large  error  magnification 
is  expected  since  the  radial  stresses  are  about  10  times  larger  than 
the  expected  stress  difference.  In  such  cases,  the  flow  equations  are 
not  sensitive  to  the  precise  magnitude  of  §,  and  $ cannot  be  accurately 
estimated  by  LASS.  Note  that  0.5  GPa  is  about  15%  of  the  peak  stress 
difference  calculated  for  LASS  1. 
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To  estimate  the  magnitude  of  errors  attributable  to  the  specific 
modeling  function  or  characterization  selected  for  the  gage  flow  surface, 
we  performed  a second  piecewise  linear  analysis  of  LASS  1 using  seven 
knots  instead  of  five  and  quadratic  exponential  rather  than  linear 
exponential  functions  for  the  particle  velocity  amplitude  link  lines. 
Because  this  analysis  introduced  small  artifacts  in  the  representation 
of  the  particle  velocity  gage  flow  surface  near  the  locus  of  local  maxima, 
this  treatment  was  somewhat  less  accurate  than  the  reported  LASS  1 
analysis.  Nevertheless,  the  differences  in  gage  flow  surface  representa- 
tions between  the  two  analyses  are  probably  typical  of  the  largest  model- 
ing variations  or  errors  possible  in  the  reported  LASS  1 analysis. 

Since  at  the  center  of  the  measurement  interval  (h  = 1.05  m)  the 
differences  in  the  strains  and  stress  difference  calculated  by  the  two 
analyses  are  about  10%  of  the  peak  values  of  the  calculated  quantities, 
this  may  be  taken  as  an  estimate  of  the  modeling  error  in  the  LASS  1 
results.  This  estimate  is  consistent  with  our  finding  reported  last 
year  of  5-10%  error  in  strains  calculated  with  piecewise  linear  LASS  on 
a test  problem  for  which  the  solution  was  known.  It  is  also  consistent 
with  the  observation  in  Section  4 that  the  crude  global  analysis  of 
LASS  1 performed  as  a feasibility  demonstration  of  the  global  technique 
gave  surprisingly  reasonable  results  in  spite  of  large  modeling  errors. 

An  estimate  of  the  total  error  in  the  calculated  results  is  then 
the  root  mean  square  of  the  experimental  and  modeling  errors.  Therefore 
the  reported  LASS  1 results  at  the  middle  of  the  h-interval  are  estimated 
to  be  accurate  to  about  ± 20%  of  their  peak  values. 


CONCLUSIONS 

Both  field  tests  yielded  a 100%  data  return  and  excellent  quality 
records.  LASS  1 data  are  slightly  superior  to  those  from  LASS  2,  as 
expected,  since  more  nearly  ideal  procedures  were  possible  in  the  first 
shot  for  logistics  reasons.  The  stress-strain  trajectories  calculated 
by  LASS  for  each  of  the  two  field  tests  are  estimated  to  be  accurate 
near  the  center  of  the  instrumented  h-interval  to  about  20%  of  their 
peak  values.  Results  from  the  two  tests  are  in  close  agreement  qualita- 
tively. Quantitatively,  also,  most  comparable  calculations  are  within 
20%  of  each  other,  indicating  agreement  between  the  two  shots  to  well 
within  experimental  error.  Efforts  were  made  both  in  advance  and  on 
site  to  avoid  differences  in  geology  between  the  two  test  regions,  which 
are  separated  by  about  60  m.  Although  such  geology  differences  have 

been  observed  at  the  Pre-Dice  Throw  II  site  in  borings  separated  by 
15 

similar  distances,  none  were  detected  or  are  known  in  the  present  case. 

Two  other  material  properties  tests  have  been  performed  in  the  Pre- 

Dice  Throw  II  wet  clay  stratum  in  which  we  fired  the  field  tests: 

* 

CIST  15  by  AFWL  and  quasi-static  laboratory  uniaxial  strain  (UX) 

16 

tests  by  WES.  The  results  of  all  three  tests  were  compared  by 
Dr.  J.  G.  Jackson,  Jr.,  of  WES  at  the  1977  DNA  strategic  structures 
meeting  and  are  summarized  here  in  a stress-volume  strain  plot  in 
Figure  41.  The  LASS  results  are  represented  by  the  equilibrium  com- 
pression curves  (see  Figure  40)  based  on  4:he  data  in  Figures  27  and  37. 

The  CIST  results  are  represented  by  a nearly  linear  mean  stress-volume 
strain  relation  fit  to  the  measured  velocity  histories  and  an  assumed  fail- 
ure surface;  and  the  UX  results  by  the  measured  uniaxial  strain  stress- 
strain  path.  All  results  pertain  to  the  wet  clay  stratum  centered  at  a 

* 

Private  communication,  G.  Wayne  Ullrich,  Air  Force  Weapons  Laboratory. 
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FIGURE  41  STRESS  VERSUS  VOLUME  STRAIN  FOR  PRE  DICE  THROW  II  WET  CLAY  FROM 
THREE  DIFFERENT  TESTS:  LASS,  CIST,  AND  UX 
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depth  of  about  3.7  m.  Because  the  strength  of  the  test  material  is 

known  a priori  to  be  negligible  on  the  stress  scale  of  Figure  41  (Ullrich 

-3  16 

in  fact  estimated  it  to  be  about  10  GPa  from  the  CIST  data),  the 

stress-strain  trajectories  for  each  of  the  three  investigations  are 

essentially  hydrostats.  However,  each  "hydrostat"  is  determined  for 

different  deformations  (spherical,  cylindrical,  and  uniaxial),  different 

4 -1  -1 

strain  rates  (from  about  10  s to  about  Is  ),  and  different  specimen 
sizes  and  environments  (laboratory  and  in  situ). 

Figure  41  shows  clearly  that  the  test  material  is  significantly 

less  compressible  (by  a factor  slightly  less  than  three)  in  the  present 

high  strain  rate  in  situ  tests  than  in  the  quasi-static  laboratory  tests. 

The  CIST  results,  determined  in  situ  and  at  intermediate  strain  rates 
2 -1 

(~10  s ),  fall  between  the  LASS  and  UX  responses.  It  should  be  noted, 
however,  that  all  the  CIST  data  are  determined  at  stresses  below  0.05 
GPa  and  that  the  accuracy  of  the  results  extrapolated  to  higher  stresses 
is  not  known. 

We  conclude  from  Figure  41  that  the  results  of  the  three  investiga- 
tions indicate  statistically  significant  response  dependence  on  loading 
path,  strain  rate,  and  environment.  Additional  sensitivity  studies 
beyond  the  scope  of  the  present  work  would  be  required  to  determine 
which  of  these  factors  is  dominant  in  altering  response.  Evidently  any 
material  model  constructed  for  the  wet  clay  test  material  may  be  in 
error  by  about  a factor  of  three  unless  it  is  based  on  data  obtained 
under  the  conditions  of  interest  or  it  contains  the  appropriate  path, 
strain  rate,  and/or  environment  dependence. 

The  objective  of  the  work  described  in  this  section  was  to  proof 
test  the  IASS  technique  using  field  data.  We  feel  that  the  proof  test 
Has  been  very  successful  and  that  the  technique  is  ready  for  use  in 
characterizing  sites  of  interest  to  DNA.  In  addition,  valuable  new 
characterization  data  have  been  provided  for  the  Pre-Dice  Throw  II  site. 

I 
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4.  ANALYSIS  DEVELOPMENTS 


LASS  has  been  successfully  proof  tested  with  field  data  and  in  the 
future  may  be  one  of  several  techniques  used  concurrently  to  determine 
material  properties  at  sites  of  interest  to  DNA.  Therefore  this  year 
we  made  extensions  to  the  technique  and  performed  validation  studies  to 
establish  and/or  increase  the  confidence  limits  for  LASS  results.  The 
purpose  is  to  provide  information  for  properly  weighting  LASS  results 
in  material  modeling  efforts.  These  extensions  and  validation  studies 
are  described  below. 

EXTENSIONS 

The  three  major  additions  to  LASS  capabilities  this  year  were; 

(1)  digital  smoothing  for  the  field  records,  (2)  an  alternative  form 
of  Lagrange  analysis  called  the  global  analysis,  and  (3)  statistical 
error  propagation  expressions  for  the  piecewise  linear,  global,  and 
finite  difference  Lagrange  analyses. 

Digital  Smoothing  of  Field  Data 

Previously  in  this  program  we  smoothed  the  field  data  "by  eye"  to 
eliminate  random  noise  and  obvious  artifacts  before  Lagrange  analysis. 

The  smoothing  was  done  by  the  experimenter  based  on  his  evaluation  of 
the  specific  data  and  on  his  past  experience  with  similar  experiments. 
Although  this  procedure  is  scientifically  valid,  it  is  not  reproducible 
or  easily  documented  for  those  who  may  wish  to  reinterpret  the  data  or 
the  LASS  results. 

This  year  we  formulated  and  reduced  to  code  a digital  low-pass 
window  filtering  technique.  The  technique  was  used  to  smooth  the  LASS  1 
and  LASS  2 field  data  (Section  3).  The  basic  approach  is  to  (1)  digitize 
each  record  at  regular  time  steps,  (2)  select  the  number  of  time  steps, 
data  points,  or  the  "window-width,"  centered  about  the  time  step  in 
question  over  which  the  data  is  to  be  smoothed,  (3)  Fourier-decompose 
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the  data  within  the  window,  (4)  filter  out  the  noise  using  an  appropriately 
shaped  filter,  (5)  retransform  to  obtain  the  amplitude  of  the  smoothed 
function  at  the  center  of  the  window,  and  (6)  repeat  at  the  next  time 
step.  The  width  of  the  window  and  the  shape  of  the  filter  are  continuously 
adjustable  as  the  window  moves.  In  this  way  we  can  account  for  the  time- 
variant  nature  of  the  signal  and  noise.  In  addition  we  constrain  the 
procedure  to  avoid  (1)  pulse-spreading  equivalent  to  information  traveling 
at  greater  than  local  sound  velocity  (causality  violation)  and  (2)  mean- 
square  amplitude  attenuation  (power  loss).  The  specific  procedures  used 
to  smooth  the  LASS  1 and  IASS  2 data  are  given  in  Appendix  C. 

Global  Lagrange  Analysis 

The  spherical  and  cylindrical  Lagrange  analyses  are  means  of 
numerically  solving  the  flow  equations  (conservation  of  mass  and  momentum) 
for  unmeasured  quantities  of  interest  starting  with  several  experimental 
measurements  of  Lagrange  radial  stress  and  particle  velocity  histories 
(see  e.g. , Appendix  A) . These  analyses  are  extremely  sensitive  to  the 
local  estimated  values  of  partial  derivatives,  which  are  determined 
from  the  discrete  spatially  separated  gage  records.  Because  these 
field  records  are  often  noisy  and  are  determined  at  widely  separated 
coordinates  in  the  flow,  the  accuracy  of  the  analysis  depends  critically 
on  the  numerical  procedures  used  to  estimate  the  required  intermediate 
quantities  from  the  individual  field  records.  Three  possible  approaches 
are  (1)  straight  finite  differences  in  which  time  derivatives  are 
taken  directly  from  digitized  gage  records  and  spatial  derivatives 
are  taken  by  comparing  adjacent  records,  (2)  the  piecewise  linear 
analysis  discussed  in  Appendix  A in  which  individual  records  are 
represented  as  piecewise  linear  functions  of  time  but  all  records 
are  used  to  establish  spatial  dependences  and  profile  amplitude  at 
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an  individual  gage  location,  and  (3)  the  global  approach  discussed 
below  in  which  a more  general  global  function  (not  usually  piecewise 
linear  in  time)  is  fit  to  all  the  data. 

With  perfect  data  a sophisticated  version  of  the  finite  difference 
approach  would  be  optimal.  However,  because  of  the  noise  in  individual 
field  records,  the  finite  difference  method  is  unsatisfactory  in  this 
application.  Therefore  in  previous  work  in  this  program  we  developed 
the  piecewise  linear  Lagrange  analysis,  which  averages  over  all  the 
gage  records.  We  have  demonstrated  here  and  in  the  earlier  work  that 
this  method  gives  good  results  for  all  calculated  quantities  if  the 
stress  range  of  the  flow  measurements  is  comparable  to  the  stress  difference 
in  the  flow.  If  not,  the  strains  can  still  be  accurately  calculated, 
but  the  flow  is  insensitive  to  the  stress  difference  which  is  not 
accurately  predicted.  To  improve  the  stress  difference  estimate  in 
this  unfavorable  condition,  we  formulated  and  reduced  to  code  the  global 
analysis,  which  makes  nearly  optimal  use  of  the  available  data. 

In  this  analysis  a global  model  function,  f , of  Lagrange  position 
♦ s 

h,  relative  time  t,  and  parameters  a^ , i = 1,2,  ...  N,  is  selected 
and  fitted  to  all  the  stress  records;  another  function,  f^,  is  fitted 
to  all  the  particle  velocity  records: 

f = f (h,t,a  ) i = 1,2,  . . . N 

s s i 

f = f (h,t,a  ) j = 1,2,  ...  M. 

v v j 

The  fitting  operation  determines  values  for  the  model  function  parameters 

a.  and  a,  to  minimize  the  Chi-square  error  between  the  model  function 

and  the  gage  records.  Once  analytic  global  model  functions  have  been  generated 


* 

Time  measured  from  the  foot  of  the  wave. 
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for  the  stress  and  particle  velocity  fields,  the  flow  equations  (Appendix 
A)  are  solvable  directly  since  the  necessary  partial  derivatives  and 
integrals  can  be  evaluated  analytically. 

Like  the  piecewise  linear  Lagrange  analysis  discussed  in  Appendix  A, 
the  global  analysis  reduces  the  effect  of  random  gage  errors  by  systemati- 
cally averaging  (globally  smoothing)  the  gage  records.  However,  it  goes 
beyond  the  piecewise  linear  analysis  in  that  it  uses  the  entire  gage 
record  rather  than  just  its  magnitude  at  characteristic  points,  i.e., 
at  the  knots  in  the  piecewise  linear  representations.  In  this  way  the 
global  analysis  makes  more  nearly  optimal  use  of  the  data  by  extracting 
information  from  the  regions  between  the  characteristic  features,  thereby 
improving  time-derivative  estimates. 

The  global  analysis  is  also  convenient  for  enforcing  known  physical 
constraints.  Such  constraints  can  be  enforced  during  global  fitting  in 
three  separate  ways: 

(1)  The  constraint  can  be  used  to  eliminate  a parameter 
from  the  global  model  before  fitting  is  attempted. 

(2)  The  constraint  can  be  strictly  enforced  during  fitting 
by  using  Lagrange  multipliers. 

(3)  The  constraint  can  be  weakly  enforced  by  adding  a penalty 
term  to  the  chi-square  sum  before  attempting  the  fit. 

Several  examples  are  discussed  below. 

Causality  is  enforced  by  insisting  that  both  f and  f be  identically 

v s 

zero  before  the  shock  front  arrives.  That  is,  the  global  model  is  designed 

so  that  f and  f are  always  zero,  when  t < 0,  where  t is  relative  time, 
v s 

(Recall  that  t = 0 on  shock  arrival.)  The  arrival  time  of  the  shock 
front  is  fitted  to  a uniform  (elastic)  or  nonuniform  (Inelastic)  form, 
as  the  case  may  be,  to  relate  relative  time  to  real  time. 

The  stress  difference  0 is  known  to  start  from  zero  at  t = 0 and 
become  positive  on  arrival  of  the  shock  front,  t > 0.  In  practice,  one 
parameter  of  the  global  model  is  eliminated  with  the  first  (t  = 0)  condition. 
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The  second  condition  (t  > 0)  is  enforced  during  fitting  with  an  inequality 
constraint . 

In  the  initial  (elastic)  regime,  the  stress  difference,  tangential 
strain,  and  radial  strain  are  linked  by  the  following  relation: 

$ = g - C > 
r t 

where  g is  constant.  When  the  experimental  data  contain  large  errors, 
it  may  prove  advantageous  to  enforce  this  constraint  and  thereby  eliminate 
another  parameter  from  the  model.  If  the  experimental  data  are  well  defined, 
the  above  condition  can  be  weakly  enforced  during  the  fitting  procedure. 

(Both  procedures  require  an  a priori  estimate  of  a constant  g.)  Alterna- 
tively, at  high  stresses  is  often  known  a priori  to  be  negligible,  and  the 

condition  <6=0  can  be  used  to  correlate  the  forms  of  f and  f . 

v s 

In  some  flows  the  peak  signal  will  occur  at  the  same  relative  time 
for  equidistant  velocity  and  stress  gages.  This  constraint  is  enforced 
by  insisting  that 


at  the  same  relative  time,  t. 

In  summary  the  global  analysis  is  an  alternative  form  of  LASS  that 
makes  nearly  optimal  use  of  the  field  data,  provides  built-in  data 
smoothing,  and  permits  enforcement  of  known  physical  constraints 
although  it  is  not  necessary  to  make  such  assumptions  on  the  material 
model.  As  described  later  in  this  section  we  have  derived  error 
propagation  analyses  for  it  and  have  used  it  in  several  validation 
studies  this  year. 
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Error  Propagation  Analyses 

This  year  we  performed  the  first  quantitative  treatment  of  error 
propagation  in  the  Lagrange  analyses.  For  the  finite  difference,  global, 
and  piecewise  linear  forms  of  LASS,  we  derived  approximate  expressions 
for  the  magnitudes  of  the  errors  in  calculated  quantities  resulting 
from  propagation  of  the  dominant  experimental  uncertainties,  namely, 
amplitude  uncertainty  in  the  gage  records  and  gage  position  uncertainties. 
The  treatment  was  statistical  with  liberal  use  of  simplifying  assumptions 
and  first-order  mathematical  approximations.  The  analyses  and  resulting 
expressions  are  presented  in  Appendixes  D through  F.  The  global  form, 
Appendix  E,  was  reduced  to  code  and  used  in  several  test  problems  as 
described  below. 

The  addition  of  error  propagation  capabilities  to  LASS  is  a significant 
advance  not  only  because  it  allows  us  to  assign  confidence  levels  to 
LASS  results,  but  also  because  it  permits  us  to  perform  pretest  sensitivity 
studies  to  optimize  field  test  design.  For  example,  proposed  tests  can 
be  evaluated  with  regard  to  both  maximum  expected  accuracy  and  cost- 
effectiveness  of  data  channel  utilization. 


VALIDATION  STUDIES 

Validation  studies  were  performed  to  refine,  determine  the 
sensitivity  of,  or  otherwise  provide  information  for  evaluating  or 
improving  the  accuracy  of  present  Lagrange  analyses.  We  studied  the 
piecewise  linear  form  of  LASS  developed  previously  and  used  to  reduce 
the  LASS  1 and  2 field  test  data,  as  well  as  the  global  form  of  LASS 
developed  this  year  and  used  to  estimate  the  errors  in  the  field  test 
results.  In  this  work  we  used  an  elastic  test  problem  for  which  we 
have  analytic  solutions  for  all  quantities  of  interest.  A more  complex 
inelastic  test  problem  for  which  we  have  numerical  solutions  for  all 
quantities  of  interest,  and  the  LASS  1 field  data.  Below  we  briefly 
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describe  these  studies  categorized  by  the  form  of  LASS  analysis  and  the 
data  source  and  give  the  primary  conclusions  drawn  from  them. 

Piecewise  Linear  IASS,  Elastic  Test  Problem 

A specific  set  of  analytic  solutions  to  an  elastic  spherical  flow 

7 

problem  was  described  by  Grady  in  a previous  report  in  this  program. 

Five  Lagrangian  radial  stress  and  radial  particle  velocity  histories 
at  coordinate!  4,  5,  6,  7,  and  8 from  the  analytic  solution  for  this 
problem  are  reproduced  in  Figure  42  (note  that  all  quantities  are 
normalized  to  be  dimensionless) . These  histories  correspond  to  error- 
free  gage  records  and  provide  the  input  data  for  testing  Lagrange 
analyses.  Grady  analyzed  these  data  by  a modified  finite  difference 
form  of  LASS  to  obtain  the  results  shown  in  Figure  43.  The  accuracy 
is  within  1/2%.  Because  the  data  are  perfect,  finite  differences  can 
give  time  derivatives  as  accurately  as  desired.  Consequently  the  errors 
that  accrue  are  primarily  from  imperfect  spatial  derivative  estimates, 
which  must  be  based  on  just  the  five  gage  stations  available.  Grady's 
work  demonstrated  that  the  Lagrange  analysis  approach  is  feasible  given 
a reasonable  number  of  realistically  separated  and  accurate  gages.  Note 
also  that  the  stress  difference,  (t>,  is  well  defined  by  the  flow  equations 
(i.e.,  the  Lagrange  analysis)  in  this  case  in  which  the  magnitude  of  <J> 
is  comparable  to  the  magnitude  of  the  radial  stress. 

We  performed  three  piecewise  linear  Lagrange  analyses  on  this 
problem  to  test  improved  numerical  techniques  and  to  determine  how  the 
LASS  accuracy  would  be  degraded  by  the  modeling  errors  introduced  in 
replacing  the  correct  profiles  with  piecewise  linear  representations. 
Basically  the  piecewise  linear  LASS  is  itself  a finite  difference 
analysis  but  with  very  large  steps,  giving  much  cruder  time  derivative 
estimates  than  Grady's  analysis.  As  stated  earlier,  this  approach  is 
used  because  conventional  finite  difference  Lagrange  analyses  of  field 
data  are  unstable  due  to  the  noise  in  the  data. 
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FIGURE  43  FINITE  DIFFERENCE  LASS  CALCULATIONS  (SOLID  CURVES)  COMPARED  WITH 
THE  ANALYTIC  PROFILES  (DASHED  CURVES)  AT  LAGRANGE  COORDINATE 
h = 6,  ELASTIC  TEST  PROBLEM 
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In  the  first  piecewise  linear  analysis,  the  "gage  records,"  i.e., 
the  data  in  Figure  42,  were  represented  with  10  knots  or  piecewise  linear 
segments,  a relatively  large  number.  The  second  analysis  used  only  two 
knots,  the  minimum  number,  maximizing  modeling  error.  The  third  again 
used  10  knots,  but  the  presumed  gage  locations  were  altered  randomly 
to  determine  analysis  sensitivity  to  the  surveying  problem  encountered 
in  the  field.  Gage  location  errors  are  expected  to  significantly  impare 
accurate  estimates  of  spatial  derivatives.  The  new  perturbed  locations 
of  the  gages  are:  3.886,  4.962,  6.037,  6.977,  and  8.063  for  particle 
velocity  and  3.899,  4.959,  6.031,  7.044,  and  7.991  for  stress. 

The  results  of  the  three  calculations,  namely,  the  calculated  strains 
and  stress  difference  at  the  middle  gage  position,  h = 6,  are  shown  in 
Figure  44  as  solid  lines.  The  analytic  solutions  are  shown  dashed.  The 
results  indicate  that  the  modeling  errors  detectably  affect  the  calculations, 
but  that  the  quantitative  and  qualitative  nature  of  the  flow  is  reproduced 
with  excellent  accuracy  even  for  the  severe  two-knot  case.  The  conclusion 
is  that  the  present  piecewise  linear  analysis  is  accurate  and  stable  with 
respect  to  modeling  and  gage  position  errors.  In  particular,  the  stress 
difference  estimates  are  quite  acceptable,  within  10%  for  10  knots  and 
accurate  records  and  positions.  This  again  is  because  4>  is  well  defined 
by  the  flow  parameters  in  regions  where  it  is  comparable  in  magnitude 
to  the  radial  stress,  at  these  attenuation  rates. 

All  three  of  these  calculations  include  two  significant  improvements 
over  previous  piecewise  linear  analyses.  First,  the  knots  in  the  stress 
and  particle  velocity  records  are  chosen  to  occur  at  the  same  relative 
times  at  a given  Lagrange  coordinate.  This  better  correlates  the  stress 
and  particle  velocity  record  representations  and  thereby  improves  the 
accuracy  of  calculations  comparing  derivatives  of  different  gage  types. 
Second,  we  determined  that  the  analysis  is  relatively  more  sensitive  to  the 
time-distance  representation  of  the  foot  of  the  wave  than  to  the  other 
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link-lines.  Therefore,  greater  attention  was  devoted  to  selection  of 
the  functional  form  for  this  link  line. 

Piecewise  Linear  LASS,  Inelastic  Test  Problem 

An  inelastic  test  problem  involving  a complex  flow  with  an  evolving 

1 

precursor  was  described  last  year.  This  problem  was  generated  with  the 
spherical  PUFF  code  using  a known  material  model.  We  showed  in  that 
work  that  piecewise  linear  LASS  could  calculate  radial  stress-volume 
strain  trajectories  accurate  to  5-10%  at  the  center  of  the  gaged  region 
in  such  a flow. 

This  year  we  reanalyzed  that  problem  using  the  two  improvements  in 
the  piecewise  linear  analysis  described  above.  The  original  and  present 
radial  stress-volume  strain  calculations  and  the  correct  trajectories 
from  spherical  PUF'F  are  shown  in  Figure  45.  In  the  new  calculation, 
the  accuracy  has  been  somewhat  improved  in  the  center  of  the  gaged 
region  and  significantly  improved  at  the  edges.  The  calculated  <t> 
histories  were  noisy,  because  <J>  is  small  compared  with  the  radial 
stresses.  On  the  basis  of  this  work  we  now  estimate  that,  with  five 
gage  stations  at  which  accurate  stress  and  particle  velocity  histories 
are  obtained,  strain  calculations  to  within  5%  are  possible  with  the 
piecewise  linear  LASS. 

Global  LASS,  Elastic  Test  Problem 

The  elastic  test  problem  was  used  to  perform  four  sensitivity 
studies  of  the  global  Lagrange  analysis.  The  intent  was  to  determine 
the  effects  on  this  analysis  of  modeling  errors  (case  1),  modeling 
errors  and  amplitude  errors  in  the  gage  records  (case  2),  modeling 
errors,  amplitude  errors,  and  gage  position  errors  (case  3),  and  fewer 
gage  stations  (three  instead  of  five)  with  the  same  errors  as  in  case  3 
(case  4).  The  results  are  shown  in  Figure  46  along  with  the  analytic 
solution  (dashed). 
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BY  PUFF  CODE 

Note  that  this  year’s  calculations  are  significantly  closer  to  the  correct  (PUFF)  trajectories 
everywhere,  but  especially  near  the  ends  of  the  interval  (h  - 0.076  m and  h = 0.113  ml 
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SOLUTIONS  (DASHED  CURVES)  AT  LAGRANGE  COORDINATE  h = 6,  ELASTIC 
TEST  PROBLEM. 

Case  1,  modeling  errors  only;  Case  2,  Case  1 plus  amplitude  errors;  Case  3,  Case  2 plus 
gage  position  errors;  Case  4,  Case  3,  but  only  the  middle  three  gages  used  in  global 
Lagrange  analysis. 
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FIGURE  46  GLOBAL  LASS  CALCULATIONS  (SOLID  CURVES)  COMPARED  WITH  ANALYTIC 
SOLUTIONS  (DASHED  CURVES)  AT  LAGRANGE  COORDINATE  h = 6,  ELASTIC 
TEST  PROBLEM  (CONCLUDED) 

Case  I,  modeling  errors  only;  Case  2,  Case  1 plus  amplitude  errors;  Case  3,  Case  2 plus 
gage  position  errors;  Case  4,  Case  3,  but  only  the  middle  three  gages  used  in  global 
Lagrange  analysis. 


Case  1 — The  stress  and  particle  velocity  gage  profiles  in  Figure  42 


were  fit  to  the  modeling  function 

2 3 

(At  + Bt  + Ct  ) exp  (Dt)/(h  + 1) 

where  t is  relative  time  at  the  Lagrange  coordinate  h.  The  fit  is 
reasonably  good  as  can  be  seen  from  the  top  two  plots  for  case  1 in 
Figure  46.  The  strains  and  stress  difference  calculated  at  h = 6 are 
shown  in  the  bottom  two  plots.  The  differences  between  calculations 
and  analytic  solutions  are  due  solely  to  modeling  errors.  Better  model 
functions  would  reduce  these  errors,  but  the  present  fit  would  be  more 
than  adequate  for  field  work. 

-4 

Case  2 — Random  amplitude  noise  with  standard  deviation  of  10  was 

added  to  the  gage  records  of  Figure  42  and  global  LASS  rerun.  The 

perturbed  records  and  new  model  functions  are  shown  in  the  first  two 

plots,  case  2,  Figure  46,  and  the  calculations  are  compared  with  the 

analytic  solutions  in  the  las  + two  plots.  In  addition,  the  propagated 

-4 

error  due  to  amplitude  uncertainties  of  10  was  calculated  by  the 
techniques  described  in  Appendix  I)  and  are  shown  in  Figure  46.  For  the 
strains,  these  are  too  small  to  be  in  this  plot. 

Case  3 — Gage  position  perturbations  described  in  the  discussion  of 
the  piecewise  linear  LASS  calculation  for  the  elastic  test  problem  were 
then  added  and  global  LASS  was  rerun  and  errors  were  reestimated.  The 
results,  the  analytic  solutions,  and  the  estimated  errors  are  s town  in 
case  3,  Figure  46.  The  differences  between  the  calculations  and  the 
analytic  solution  are  due  to  modeling,  position,  and  gage  location  errors. 
The  error  estimate  does  not  include  modeling  errors. 

Case  4 — The  curves  presented  in  Cases  1 through  3 were  obtained  by 
fitting  simulated  data  from  five  particle  velocity  and  five  stress  gages. 
When  only  three  particle  velocity  and  stress  gages  were  used,  the  results 
presented  in  case  4 were  obtained.  Although  profiles  are  shown  at  all  five 
gage  stations  in  the  first  two  plots  of  case  4,  only  the  middle  three  were 
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used  in  the  global  analysis.  Note  that  the  assigned  errors  increase  by 
a factor  of  nearly  2 when  fewer  gage  records  are  available  for  global 
fitting. 

These  calculations  demonstrate  the  type  of  sensitivity  studies 
that  can  be  performed  using  a LASS  analysis  with  error  estimation 
capabilities.  For  example,  before  a field  test,  we  could  specify  the 
degree  of  the  maximum  allowable  uncertainty  acceptable  in  the  results 
and  the  expected  precision  of  the  experimental  measurements  and  then 
determine  how  many  data  channels  were  required  to  achieve  the  specified 
accuracy.  We  could  examine  the  cost  effectiveness  of  deleting  channels, 
improving  gage  surveying  procedures,  and  so  on.  The  present  results 
also  demonstrate  that  the  global  LASS  is  highly  stable  to  noise  and 
position  errors. 

Global  LASS,  LASS  1 Field  Test 

A preliminary  global  analysis  of  the  LASS  1 field  data  presented  in 
Section  2 was  performed  to  test  this  analysis  with  real  data  and  to 
generate  approximate  error  assignments  for  the  Lagrange  analysis  results 
in  Section  3.  Although  a very  crude  model  producing  severe  modeling 
errors  was  used  in  this  preliminary  analysis,  plausible  estimates  of 
strains  and  stress  difference  were  obtained  and  noise  and  gage  position 
errors  were  properly  assigned. 

The  stress  and  particle  velocity  profiles  were  fit  to  the  function 

2 2 
At  exp  (Bt  + Ct  + Dh) 

where  t is  relative  time  at  the  Lagrange  coordinate  h.  This  form  was 
selected  for  the  feasibility  demonstration  for  convenience  and  a semblance 
of  plausibility.  Other  forms  can  be  accommodated  with  minor  changes; 
however,  the  noise  and  gage  position  errors  are  propagated  correctly 
regardless  of  the  selected  form.  As  discussed  in  Section  3,  we  made 
the  following  error  assignments: 


r 


particle  velocity  ± 5 ms 
radial  stress  ± 0.013  GPa 
gage  position  ± 38  mm 

The  results,  shown  in  Figure  47,  are  quite  similar  to  the  piecewise 
linear  analysis  presented  in  Section  3 in  spite  of  the  large  modeling 
errors  evident  in  plots  (a).  The  most  important  results  in  this  figure 
are  the  ratios  of  estimated  errors  to  the  calculated  quantities.  These 
relative  error  estimates  were  used  in  the  LASS  analysis  described  in 
Section  3. 

Finite  Difference  LASS,  Elastic  Test  Problem  and  LASS  1 

Modified  finite  difference  Lagrange  analyses  of  both  the  elastic 

problem  and  the  LASS  1 field  data  were  performed  to  determine  the 

stability  of  an  analysis  using  a nonglobal  approach  to  evaluating  the 

temporal  derivatives  in  the  flow  equations.  The  spatial  derivatives  were 

still  determined  by  path  lines  and  hence  use  global  smoothing.  This 

modified  finite  difference  approach  is  the  method  used  by  Grady  to 

obtain  the  elastic  results  shown  in  Figure  43.  The  elastic  profiles 

with  both  noise  and  gage  location  errors  (case  3,  Figure  46)  were  used. 

The  results  for  both  the  elastic  problem  and  the  LASS  1 field  data  are 

presented  in  Figure  48.  The  radial  strains  are  smooth  (since  they 

involve  only  integrals  and  not  derivatives  of  the  gage  records)  while 

the  volume  and  radial  strains  are  slightly  unstable.  The  stress 

difference,  appearing  in  the  flow  equations  as  a small  difference  in  two 

large  derivatives  of  the  experimental  data,  is  extremely  unstable  and 

oscillates  randomly  off-scale  for  the  elastic  problem.  The  scale  on 

which  the  LASS  calculation  is  shown  is  very  insensitive  compared  with 

-3 

the  expected  stress  difference  amplitude  of  ~10  GPa  or  it  would  look 
the  same.  These  results  demonstrate  the  significant  improvements 
obtained  by  the  global  and  piecewise  linear  Lagrange  analysis  approaches 
when  the  experimental  data  contain  errors. 
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(a)  ROUGH  FITS  OF  SOME  OF  THE  GLOBAL  MODEL  FUNCTIONS  (DASHED)  TO 
SOME  OF  THE  FIELD  TEST  DATA. 


(b)  GLOBAL  LASS  CALCULATIONS  OF  STRAINS  AND  STRESS  DIFFERENCE  WITH 
ERROR  ESTIMATES,  FOR  A PARTICLE  ENTERING  THE  FLOW  AT  ABOUT  2.5  ms. 


(c)  RADIAL  STRESS-VOLUME  STRAIN 

AT  THE  SAME  PARTICLE. 


MA-4663-50 


FIGURE  47  PRELIMINARY  GLOBAL  ANALYSIS  AND  ERROR  ESTIMATES  BASED  ON  LASS  1 

field  test  data.  ° 

These  are  preliminary  results  intended  to  test  and  demonstrate  the  global  LASS  analysis 
and  should  not  be  quoted  or  used  in  any  quantitative  research  related  to  LASS  1 and 
LASS  2 field  tests 
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(a)  ELASTIC  TEST  PROBLEM  WITH  NOISE. 


(b)  LASS  1 FIELD  DATA 
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FIGURE  48  FINITE  DIFFERENCE  LASS  CALCULATIONS  OF  ELASTIC  TEST  PROBLEM  WITH 

NOISE  AND  OF  LASS  1 FIELD  DATA  FOR  COMPARISON  WITH  PIECEWISE  LINEAR 
AND  GLOBAL  CALCULATIONS. 

Note  that  the  stress  difference  calculations  are  very  noisy  (off  scale  for  case  a)  and  are 
reported  at  smaller  Lagrange  coordinates  than  the  strain  calculations. 
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CONCLUSIONS 


The  most  important  conclusions  of  the  validation  studies  were  that 
(1)  stress  differences  can  be  calculated  by  current  techniques  to  10% 
given  accurate  gage  records  obtained  in  a flow  region  where  the  radial 
stress  is  comparable  to  the  stress  difference,  (2)  strains  can  be 
calculated  to  within  5%  given  a realistic  number  of  accurate  records 
even  from  a complex  flow,  (3)  the  global  and  piecewise  linear  Lagrange 
analyses  are  not  unstable  to  reasonable  perturbations  (errors)  in  the 
experimental  data,  and  (4)  the  new  error  propagation  capabilities  can 
be  used  to  assign  confidence  levels  to  LASS  results  or  to  perform 
sensitivity  studies  in  advance  of  field  tests  to  evaluate  or  improve 
the  probability  of  achieving  the  desired  test  objective. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  year  LASS  has  been  successfully  proof  tested  in  conjunction 
with  two  spherical  HE  field  tests  we  performed  at  the  Pre-Dice  Throw  II 
site.  In  addition,  error  propagation  capabilities  have  been  added  to 
the  technique  and  tested.  We  estimate  that  the  LASS  results  reported 
for  the  two  field  tests  are  accurate  to  about  20%  of  the  peak  values 
of  the  calculated  quantities.  We  have  also  demonstrated  that  the  LASS 
approach  now  allows  accurate  determination  of  stress  difference  if 
reasonable  experimental  measurements  are  obtained  in  a stress  range 
comparable  to  the  stress  difference  supportable  by  the  test  material. 

Based  on  these  results  and  our  understanding  of  the  major  material 
modeling  problems  of  concern  to  DNA  and  the  status  of  techniques  available 
in  the  community  to  solve  these  problems,  we  suggest  four  research  tasks 
warranting  consideration  for  future  work. 


(1)  An  experimental  program  to  improve  dynamic  stress  measurement 
capabilities  in  the  important  stress  range  between  10“ 3 and 
10”!  GPa.  The  strengths  of  many  geologic  materials  of  interest 
lie  within  this  range  where  we  are  currently  unable  to  make 
reliable  stress  measurements.  If  reliable  measurements  are 
obtained,  LASS  and  perhaps  other  techniques  could  be  used 

to  obtain  material  properties  data  in  this  important  and 
poorly  understood  region. 

(2)  The  application  of  LASS  techniques  to  future  CISTs.  If 
the  stress  measurement  problems  can  be  solved,  LASS  can  be 
used  to  reduce  CIST  data  to  determine  unique  stress-strain 
trajectories  and  assign  confidence  levels  to  them.  This 
would  support  the  objectives  and  enhance  the  usefulness  of 
future  CISTs. 

(3)  LASS  sensitivity  analyses  of  proposed  field  tests.  Such 
analyses  can  be  used  to  establish  the  accuracy  expected  from 

a given  test  for  comparison  with  the  test  objectives.  In  this 


way  the  test  can  be  designed  either  for  optimum  accuracy 
given  testing  constraints  or  for  maximum  cost-effectiveness 
within  predetermined  accuracy  specifications. 

(4)  A sensitivity  study  to  evaluate  the  error  caused  by  applying 
LASS  to  data  from  anisotropic  or  nonuniform  sites  that  do 
not  support  strictly  one-dimensional  flows.  Since  LASS  has 
been  successfully  proof  tested,  it  is  desirable  to  remove, 
reduce,  or  quantify  this  one  restriction  on  its  applicability. 
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Appendix  A 

LAGRANGE  ANALYSIS  FOR  SPHERICAL  FLOWS 

The  LASS  analysis  for  one-dimensional  divergent  flows  calculates 
the  stress  and  strain  trajectories  of  dynamically  loaded  in  situ  site 
material  from  field  measurements  of  radial  stress  and  radial  particle 
velocity.  The  approach  is  to  use  the  field  measurements  to  numerically 
solve  the  partial  differential  equations  defining  strains  and  expressing 
the  conservation  of  mass  and  momentum  for  the  flow.  The  importance  of 
the  LASS  method  is  that  it  is  mathematically  complete,  so  no  assumptions 
whatsoever  are  required  concerning  the  thermodynamic  or  mechanical  consti- 
tutive relations  of  the  material.  Accuracy  is  limited  only  by  the  degree 
to  which  the  flow  satisfies  the  symmetry  conditions,  by  the  validity  of 
the  numerical  estimates  of  partial  derivatives,  and  by  the  experimental 
errors  in  the  field  measurements.  Thus  in  situ  HE  events  that  induce 
high  symmetry  (planar,  spherical,  or  cylindrical)  flows  serve  as  large 
material  loading  machines,  and  the  LASS  analysis  fully  defines  the  resulting 
load-response  behavior. 

The  partial  differential  equations  expressing  conservation  of  mass 

* 

and  momentum  for  planar,  spherical,  and  cylindrical  one-dimensional 
flows  may  be  written 


mass  (A-l) 


momentum  (A-2) 


♦ 

Flows  parametrized  by  one  spatial  coordinate. 


A-l 


where: 


stresses  are  taken  as  positive  in  compression 

D and  0 are  the  initial  and  present  densities 
o 

r is  the  Eulerian  radial  coordinate 

h is  the  Lagrange  coordinate  of  a particle,  i.e,,  its 
initial  radial  coordinate 

N is  a constant  with  value  0,  1,  or  2 for  planar,  cylindrical, 
or  spherical  flows,  respectively. 

t is  time 


u and  tT  are  radial  particle  velocity  and  radial  stress 
r r 

$ is  the  difference  between  radial  and  tangential  stress  (O’  - O ) 

r t 

The  definition  of  radial  particle  velocity  as  the  time  rate  of  change 
of  Eulerian  particle  position  provides  a third  flow  equation,  which 
may  be  written  in  differential  or  integral  form 


u = 
r 


( — \ 

1 dt  / 

\ / h 


(A-3) 


r = h + 


To(h) 


u (h.T)  dT 
r 


(A-3b) 


Where  T^(h)  is  the  time  of  first  motion  of  the  particle  h. 


For  N = 0 (planar  flow),  the  system  of  three  equations  and  four  flow 

variables  CT  , u , 0,  and  r can  be  solved  numerically  if  either  stress 
r r 

or  particle  velocity  histories  along  several  particle  paths  are  measured. 
For  the  divergent  flow  cases  of  interest  in  this  program  where  N = 1 or 
2 (cylindrical  or  spherical  flows),  the  momentum  equation  contains  an 
extra  term  involving  the  additional  flow  variable,  stress  difference,  $. 
For  these  flows,  both  stress  and  particle  velocity  histories  along  several 
particle  paths  must  be  measured  to  allow  the  system  of  equations  to  be 
solved  for  the  remaining  unmeasured  flow  variables. 


Note  that  for  the  divergent  flows  not  all  of  the  five  flow  variables 


are  independent.  Consequently,  measurement  of  some  pairs  of  variables 

is  not  sufficient  to  solve  for  the  others.  Inspection  of  the  equations 

shows  that  any  of  the  seven  pairs  that  can  be  formed  from  (T  , $ , and 

r 

one  of  the  remaining  displacement-related  variables  is  sufficient, 
although  some  are  computationally  or  practically  more  convenient.  Note, 
in  particular,  that  no  pair  of  displacement-related  variables  is  inde- 
pendent. All  work  to  date  has  used  the  pair  Cr  and  u . 

r r 

To  provide  stress-strain  trajectories  for  the  divergent  flows,  the 

analysis  further  calculates  strain  histories  from  the  displacement  field 

r(h,t).  The  following  principal  engineering  strain  definitions  are  used 

to  evaluate  radial  strain  € , tangential  (hoop)  strain  € , and  volume 

r t 

strain  C : 

v 


€ = 1 
r 


C = 
t 


1 


r 

h 


e = 1 

V 


= 1 


(A-4) 


(A-5) 


(A-6) 


LASS  methods  for  solving  equation  (A-l)  through  (A-6)  have  been 

1-3,7 

described  in  several  of  our  reports.  The  general  procedure  is  that 

the  field  data  are  smoothed  by  low-pass  filtering,  if  necessary,  and 
fitted  to  a functional  form.  The  functions  may  be  either  piecewise- 
linear  in  time  (the  piecewise  linear  analysis)  or  a general  analytic 
function  oi  h and  t (the  global  analysis).  Alternatively,  one 
may  estimate  the  required  partial  derivatives  by  straight  finite  differ- 
ences (the  finite  difference  analysis);  this  analysis,  however,  is 
generally  not  useful  with  field  data  because  unavoidable  experimental 
errors  in  the  data  make  the  derivative  estimates  excessively  noisy.  The 


A-3 


fitted  representations  are  then  used  to  solve  the  flow  equations,  (A-l)  - 
(A-3),  and  the  defining  strain  equations,  (A-4)  - (A-6),  for  the  instan- 
taneous values  of  all  components  of  the  stress  and  strain  tensors  at  any 
location  in  the  flow,  (Note,  however,  that  accuracy  is  greatest  near 
the  center  of  the  measurement  interval.)  Uncertainties  may  be  estimated 
by  normal  error  propagation  methods  (see  Section  4)  once  the  uncertainties 
in  the  initial  measurements  are  specified. 

The  basic  approach  of  the  piecewise  linear  LASS  analysis  developed 
in  this  program  is  described  below.  The  global  and  finite  difference 
analyses  use  standard  methods  and  are  not  discussed  here. 

The  first  step  in  the  piecewise  linear  analysis  is  to  use  the 

individual  field  records  (T  (h  , t)  and  u (h  ,t)  to  create  analytic 

r i r i 

representations  for  the  gage  flow  surfaces,  the  CT  and  u surfaces 

r r 

over  the  h-t  domain  spanned  by  the  measurements.  The  field  records 
are  constant  -h  traces  on  these  surfaces  as  shown  schematically  in 
Figure  A-l. 

The  major  problem  in  creating  accurate  representations  of  the  gage 
flow  surfaces,  or  in  performing  any  IASS  analysis,  is  the  high  level 
of  noise  and  other  anomalies  often  present  in  field  records.  These  occur 
both  because  of  the  difficulty  in  making  precise  measurements  outside 
the  laboratory  and  because  of  the  numerous  geologic  inhomogeneities 
present  in  real  sites.  To  handle  the  problem  of  errors  in  individual 
gage  measurements,  we  use  all  the  gage  measurements  to  establish  an 
averaged  flow  surface  amplitude  at  specific  h^,t  coordinates.  Having 
created  an  averaged  mathematical  representation  of  the  gage  flow  surface, 
we  can  then  generate  smoothed  estimates  of  individual  records  as  needed 
to  solve  equations  (A-l)  through  (A-6).  This  technique  is  more  circuitous 
than  using  the  field  records  directly,  but  it  significantly  improves  the 
accuracy  of  the  calculations  and  is  computationally  convenient. 
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FIGURE  A 2 


PIECEWISE  LINEAR  REPRESENTATION  OF  GAGE  RECORD. 
Points  shown  are  parameters  that  were  selected  as  best  characterizing 
the  data. 
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FIGURE  A 1 


GAGE  FLOW  SURFACE. 

The  traces  on  this  surface  provided  by  the  gage  profiles  are  used  to 
approximate  the  surface. 
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The  procedure  for  creating  the  two  averaged  gage  flow  surfaces 
consists  of  three  steps.  First,  the  significant  flow  features  on  each 
record  are  identified  and  labeled.  Significant  features  are  those 
clearly  identifiable  in  the  majority  of  the  records,  such  as  precursor 
foot,  a plastic  wave  peak,  a point  at  which  decay  rate  changes  rapidly, 
or  simply  a specific  percentage  of  peak  amplitude  for  a smooth  featureless 
profile.  Second,  since  other  features  present  in  the  individual  profiles 
are  not  deemed  statistically  significant,  the  labelled  features  on  each 
profile  are  connected  by  straight  line  segments  to  produce  piecewise 
linear  (in  time)  gage  record  representations  as  shown  in  Figure  A-2. 

The  third  step  provides  interpolation  between  Lagrange  measurements 
locations  h.  and  averages  the  individual  profiles.  This  is  accomplished 
by  linking  the  amplitude  coordinates  and  the  time  coordinates  at  each 
discrete  position  h.  for  each  feature,  such  as  peak  stress,  with  smooth 
functions  of  h.  A hypothetical  peak  stress  amplitude  versus  h linking 
is  shown  in  Figure  A-3.  Of  course,  a time  versus  h linking  for  this 
feature  is  also  necessary  to  complete  specification  of  the  peak  stress 
trace  on  the  gage  flow  surface.  Linear  or  quadratic  exponential  functions 
are  generally  used  for  the  h-linking: 


o+3h 

e , linear  exponential  fit 

yt6h+Xh  , quadratic  exponential  fit 
e 

where  O'  and  3 are  constants  created  by  the  linear  least  squares  fit, 
and  V,  6,  and  A.  are  constants  created  by  the  quadratic  least  squares 
fit.  A different  set  (0,3)  or  (y,6,A)  is  generated  for  each  amplitude 
and  each  time  link  line. 


amplitude 

or 

time 


The  fitted  link  lines  in  conjunction  with  the  piecewise  linear  form 
for  gage  records  generate  the  gage  flow  surface  representations  shown 
schematically  in  Figure  A-4.  From  this  point  solution  of  equations 
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FIGURE  A-3  AMPLITUDES  LINK  LINE  FOR  PEAK  STRESS. 

This  is  an  example  of  a link  line  generated  to  describe  the 
spatial  variation  of  one  of  the  characteristic  points  (in  this 
case,  peak  stress)  from  the  piecewise  linear  representation 
of  the  stress  records. 


FIGURE  A 4 GAGE  FLOW  SURFACE  REPRESENTATION. 

Determined  from  the  experimental  data  by  selecting  parameters  that  best  describe  the 
character  of  the  flow  within  the  quality  of  the  data.  Note  that  gage  records  are  piecewise 
linear  and  that  link  lines  are  smooth. 
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Appendix  B 
BACKGROUND 

Previous  work  by  SRI  in  the  in  situ  properties  program  and  in  re- 
lated areas  has  been  described  in  a series  of  DNA  reports  and  is  sum- 
marized  here.  In  1971  Keough  et  al.  investigated  the  feasibility  of 
applying  the  Lagrange  analysis  concept  to  radial  flows  and  of  making 
the  in  situ  ground  motion  measurements  required  for  the  analysis.  They 
presented  the  spherical  flow  equations  and  formulated  an  analysis  for 
calculating  specific  volume  from  multiple  Lagrange  histories  of  radial 
particle  velocity.  This  analysis  was  tested  against  an  analytic  solu- 
tion for  a specific  decaying  spherical  shock  problem  and  shown  to  be 
accurate.  They  also  showed  that  the  ytterbium  piezoresistive  stress 
gage  could  be  modified  for  field  use  in  divergent  flows,  and  they  sur- 
veyed possible  radial  particle  velocity  measurement  techniques. 

0 

In  1972  Smith  et  al.  modified  the  ytterbium  stress  gage  design 
and  emplacement  technique  and  obtained  successful  field  measurements  of 
stress  in  Middle  Gust,  Event  1.  The  mutual  inductance  particle  velocity 
(MIPV)  gage  concept  originally  proposed  by  Engineering  Physics  Company 
was  adopted  following  a survey  of  particle  velocity  measurement 
techniques.  Two  forms  of  this  gage  were  built,  evaluated,  and  shown  to 
be  feasible. 

7 

In  1973  Grady  et  al . developed  the  technology  and  supporting 
theory  for  the  longitudinal  MIPV  gage  currently  used  in  the  LASS  tech- 
nique and  performed  successful  field  tests  with  it.  They  further  im- 
proved ytterbium  stress  gage  response,  especially  stretch  sensitivity, 
by  modifying  the  package  design  and  emplacement  procedures.  High  ex- 
plosive proof  tests  of  stress  and  particle  velocity  gages  were  conducted 
under  conditions  similar  to  field  events,  and  stress  measurements  were 
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made  on  Middle  Gust,  Event  III.  A practical  numerical  LASS  analysis 
was  developed  and  successfully  applied  to  a spherical  elastic  problem 
with  a known  analytic  solution.  When  applied  to  experimental  records 
obtained  from  small  spherical  shots  in  Westerly  granite,  the  analysis 
gave  poorer  results,  especially  during  loading;  however,  these  results 
provided  insight  for  subsequent  analysis  development. 

O 

In  1974  Grady  et  al.  developed  the  basic  piecewise  linear  form  of 
the  numerical  LASS  analysis  now  used  for  reducing  field  quality  stress 
and  particle  velocity  records  to  dynamic  in  situ  stress-strain  trajec- 
tories. A rate- independent  variable-modulus  model  was  fitted  to  the 

17 

analysis  results  using  data  from  small  shots  in  Mixed  Company  sandstone. 

This  procedure  demonstrated  how  the  Lagrange  analysis  might  be  used  to 
support  ground  motion  calculations.  The  model,  in  conjunction  with  a 
spherical  wave  propagation  code,  produced  wave  profiles  that  compared 
favorably  with  the  experimental  data.  The  MIPV  gage  was.  developed  to 
field  instrument  capability.  Stress  and  particle  velocity  data  were 
obtained  for  the  CENSE  I event,  near  Grand  Junction,  Colorado,  and  for 
four  large  HE  tests  for  the  Essex  program  at  Fort  Polk,  Louisiana. 

In  late  1974,  Smith  et  al.**  refined  the  analysis  to  handle  more 
accurately  flows  that  exhibit  rapidly  changing  or  developing  features 
such  as  precursor  formation.  Analysis  sensitivity  to  internal  data 
handling  procedures  was  identified  during  reduction  of  records  from 
small  shots,  and  suggestions  for  improvement  were  offered.  The  stress 
and  particle  velocity  records  obtained  in  the  Essex  I,  Phase  2,  event 
in  Fort  Polk,  Louisiana,  added  to  our  experience  with  ground  motion 
measurements . 

In  1976  Rosenberg  et  al.'*’  reported  the  first  application  of  LASS 
techniques  to  field  data.  The  data  were  obtained  by  appending  ytterbium 
stress  gages  and  mutual  inductance  particle  velocity  gages  to  several 
fully  buried  shots  in  the  WES  CENSE  II  series  near  Hattiesburg,  Mississippi. 
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The  radial  stress-volume  strain  trajectory  calculated  at  the  center  of 
the  measurement  interval  agreed  with  quasi-static  laboratory  data  in 
strain  amplitude  at  peak  stress,  but  indicated  considerably  stiffer 
local  response  (higher  tangent  modulus)  for  the  dynamically  loaded  in 
situ  Material.  They  also  reported  simple  two-dimensional  code  calcula- 
tions, indicating  that  mutual  inductance  particle  velocity  gages  fielded 
inside  and  near  the  bottom  of  radial  grout  columns  would  equilibrate  to 
velocities  acceptably  close  to  the  free  field  velocity  even  with  slip 
and  significant  impedance  mismatches  between  the  grout  and  soil. 

Rosenberg  et  al.'*'  also  examined  the  accuracy  of  numerical  n.'i'fhods 
and  especially  finite  difference  approximations  used  in  the  LASS  analysis 
and  reported  that,  given  accurate  records  from  5 to  8 realistically 
separated  gage  stations,  strains  could  be  calculated  to  within  5%  to  10% 
for  complex  flows.  They  concluded  that  further  work  was  necessary  to 
calculate  stress  differences  using  LASS  methods  but  that  LASS  was  ready 
for  field  use  at  appropriate  sites. 

At  DNA's  request  SRI  has  fielded  the  gages  developed  for  the  in 
situ  program,  particularly  the  ytterbium  stress  gage,  in  several  other 
programs  during  the  past  few  years.  In  1971  Smith  described  measurements 

g 

on  two  1000-pound  high  explosive  shots  at  the  Nevada  Test  Site  (NTS). 

9 

Subsequent  ytterbium  calibration  work  by  Ginsberg  et  al.  led  Smith  to 
replot  some  of  the  waveforms:  these  are  included  in  a summary  report 
by  Jones  and  Green.'*'0  In  1972  Smith  reported  measurements  on  another 
1000-pound  test  at  NTS,  the  Mine  Dust  high  explosive  test.1^  Predictions 
and  measurements  were  in  generally  good  agreement  in  all  of  these  high 
explosive  tests  and  helped  to  establish  the  importance  of  air-filled 
porosity  in  controlling  the  rate  of  stress  wave  attenuation  in  tuff. 

In  1973  Smith  presented  ytterbium  gage  measurements  made  on  two 

12 

nuclear  shots,  Dido  Queen  and  Husky  Ace.  The  records  were  very  clean 
and  consistent,  the  best  that  had  been  obtained  at  the  time. 


I 
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In  1975  Smith!-*  reported  ytterbium  gage  measurements  made  at  the 
AFWL  Hardpan  test  site  in  a series  of  shots  designed  to  investigate  the 
stress  and  motion  fields  generated  under  multiple  loading  by  arrays  of 
HE  charges.  The  high  quality,  long  duration  records  obtained  clearly 
showed  signals  from  individual  charges.  The  measurements  demonstrated 
that  the  first  wave  produced  unexpectedly  large  alterations  in  the 
medium,  resulting  in  significant  enhancement  of  the  amplitude  of  sub- 
sequent waves. 

14 

I In  1976  Murri  and  Smith  obtained  high  quality,  long  duration 

I 

stress  measurements  in  the  0. 1-1.0  GPa  range  in  the  100-ton  and  120- ton 

if 

(90.7  and  109  Mg)  Pre-Dice  Throw  Events.  These  measurements  provided 
information  for  use  in  designing  the  600-ton  (544  Mg)  Dice  Throw  event 
ard  made  stress- time  and  stress-range  data  for  the  Pre-Dice  Throw  site 
material  available  to  the  ground  shock  community. 


I 
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Appendix  C 

DIGITAL  SMOOTHING  OF  FIELD  DATA 


D.  G.  Falconer 

INTRODUCTION 

Field  measurements  of  stress  and  particle  velocity  histories  in 
explosively  loaded  geologic  materials  contain  experimental  noise  that 
complicates  the  interpretation  and  use  of  such  data  in  Lagrange  stress- 
strain  analyses.  To  suppress  such  noise  and  thus  improve  the  scientific 
value  of  the  field  experiments,  one  can  smooth  the  gathered  data  by 
some  reproducible  and  objective  means.  The  simplest  smoothing  procedure 
is  window  averaging.  Here  a prescribed  period,  T,  is  established  and 
a simple  average  formed  about  each  data  point: 

+T/2 

v(t  - t ) dt7 

-T/2 

where  v(t)  represents  the  history  being  smoothed,  e.g.,  a particle 
velocity  record.  The  selected  period  T should  be  long  enough  to 
average  the  noise  away,  but  short  enough  to  preserve  important  physical 
features.  A sophisticated  form  of  this  approach  called  low-pass  window 
filtering  was  used  in  the  present  program. 

Spline  functions  are  another  method  of  data  smoothing  that  might  be 
applied.  For  example,  spline  functions  can  be  developed  for  the  velocity 
and  stress  records,  with  the  "knots"  located  at  the  discontinuities  in 
the  shock-wave  histories.  The  advantages  of  spline  functions  are  that 
they:  (1)  allow  for  some  curvature  in  the  shock-wave  histories  between 

knots,  (2)  provide  continuous  derivatives  for  calculating  accelerations 
and  gradients,  and  (3)  make  no  pretense  about  an  underlying  model. 


C-l 


(Piecewi se-linear  approximations  to  the  velocity  and  stress  histories  may 
be  interpreted  as  a first -order  spline  analysis.)  Two-dimensional  spline 
functions  could  provide  velocity  and  stress  estimates  over  h-t  space, 
thus  gaining  greater  averaging  power  and  eliminating  the  need  for  path- 
line approximations  for  gradient  derivatives.  This  approach  was  not 
considered  to  be  cost-effective  in  the  present  application. 

Global  modeling  is  another  approach  to  the  data-smoothing  problem. 
This  procedure  features  the  same  advantages  as  the  spline  functions 
but  allows  greater  use  of  a priori  knowledge  and  physical  law.  For 
example,  it  is  known  a priori  that  the  stress  difference  becomes  positive 
on  arrival  of  the  shock  front.  Similarly,  physical  law  requires  that  the 
shock  histories  obey  causality.  Like  the  spline  analyses,  the  global 
model  of  the  shock  histories  must  have  "knots''  between  physically  distinct 
regions  such  as  compression  and  release.  In  addition,  the  selected 
model  must  have  two  modes:  one  for  the  velocity  data  and  one  for  the 
stress  data.  Because  such  models  (1)  approach  the  data  globally,  (2) 
allow  for  a priori  knowledge  and  physical  law,  and  (3)  simplify  error 
assignments,  they  represent,  in  theory,  a best -approach  to  the  data- 
analysis  problem. 

LOW -PASS  FILTERING 

Although  window  averaging  is  quite  simple,  it  is  far  from  optimum. 

The  chief  criticism  of  this  technique  is  that  it  makes  no  provision  for 
the  time-variant  nature  of  the  signals  and  noise  in  velocity  and  stress 
records.  Also,  the  window  shape  is  uniform,  rather  than  specially 
tailored  for  the  signal  and  noise.  If  the  power  spectra  of  the  expected 
signal  and  the  interferring  noise  are  known  a priori  at  each  point  t 
in  the  time  histories,  then  an  optimum  smoothing  filter  can  be  designed. 

In  practice,  these  spectra  are  known  only  approximately,  so  that  sub- 
optimal  filtering  must  be  considered. 


A good  low-pass  filter  can  be  designed  using  the  following  steps: 

(1)  The  signal  history  is  first  assigned  a window  width  T(t)  that 
equals,  approximately,  the  period  over  which  the  signal  is 
believed  stationary — that  is,  its  power  spectrum  remains 
constant . 

(2)  The  window  shape  is  then  prescribed  by  the  Wiener  formula: 

W(f)  = A/ ( 1 + P /P  ) 
n s 

where  P (P  ) is  the  power  spectrum  of  the  signal  (noise)  at 
s n 

the  temporal  frequency  f.  Here  W(f)  is  the  Fourier  transform 
of  the  window  function  w(t), 

(3)  The  constant  A is  given  by: 


A = 1 + P (0)/P  (0) 
n s 


This  requirement  keeps  the  filtering  process  from  lowering 
the  general  signal  level. 

(4)  The  smoothed  data  v(t)  is  then  obtained  with  the  following 
formula : 

+T  ( t ) / 2 

w( t ' ) v(t  - t ' ) dt ' 

-T  ( t ) / 2 

These  steps  lead  to  a time-variant,  Wiener  filter  with  a normalized 


response. 


LASS  SMOOTHING 


A modified  version  of  the  above  procedure  was  used  to  smooth  the 
LASS  data.  Because  the  stationary  period,  T(t),  of  the  shock  histories 
was  known  only  crudely,  a simple  linear  form  was  presumed.  That  is, 


T(t)  = 


0,  t < t 


peak 


B<*  - W’  4>tpeak 


where  the  constant  B = 0.4  Us  . In  other  words,  no  filtering  was 
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carried  out  before  the  signal  peak.  After  the  signal  peak,  the  stationary 
period  (window  width)  increased  linearly  with  time. 

The  filter  form  was  presumed  to  have  a gaussian  form: 

W(f)  = exp(-f 2/f  2) 
o 

The  half-power  point  was  defined  with: 

f = 4/T 
o 

The  averaging-window  was  thus  decreased  inversely  with  the  averaging 
time,  T.  The  filtering  process  was  made  casual  by  insisting  that  the 
window  function  w(t)  be  zero  for  t < 0.  This  constraint  was  enforced 
by  multiplying  the  window  function  by  a unit  step  function  h(t).  The 
new  window  function,  w(t)h(t),  was  then  inverse  transformed  for  a new 
(and  thus  casual)  filter  function  W(f).  The  mean  and  variance  of  the 
filtered  signal  were  then  set  equal  to  the  mean  and  variance  of  the 
original  signal  to  obtain  the  final  smoothed  profiles. 
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Appendix  D 


STRESS-STRAIN  ERROR  ASSIGNMENTS  FOR  FINITE-DIFFERENCE 
LAGRANGE  ANALYSIS 

D.  G.  Falconer 


INTRODUCTION 

Field  measurements  of  ground  motion  contain  several  types  of  ex- 
perimental interferences.  Two  are  particularly  important:  amplifier 
noise  and  position  errors.  Amplifier  noise  here  refers  to  the  total  (presumably) 
additive,  zero-mean,  gaussian,  random  disturbances  in  the  recorded 
signals.  It  may  be  caused  by  site  inhomogeneities  or  many  other  factors. 
Similarly,  position  errors  refer  to  additive,  zero-mean,  gaussian, 
random  errors  in  the  geometrical  positions  of  the  velocity  and  stress 
gages.  The  principal  objective  of  this  appendix  is  to  assign  approximate 
errors  to  stress  and  strain  trajectories  calculated  by  a finite-difference 
form  of  Lagrange  analysis.  These  trajectories  are  derived  from  shock- 
wave  data  recoried  at  various  times  and  positions.  Accordingly,  the 
amplifer  noise  and  position  errors  must  be  propagated  from  the  experi- 
mental domain  to  the  stress-strain  domain. 

Error  assignments  are  made  assuming  a finite-difference  analysis 
in  which  partial  derivatives  are  calculated  along  path  lines  connecting 
corresponding  features  in  h-t  space.  (Here  h represents  the  initial 
gage  position  and  t the  elapsed  time.)  Only  first-order  approximations 
for  the  error  assignments  are  attempted.  Thus,  minor  sources  of  error 
are  omitted  and  mathematical  approximations  used  liberally.  These  con- 
ditions greatly  simplify  the  error  analysis  yet  give  satisfactory  error 
assignments  on  stress  and  strain. 

Data  smoothing  may  be  expected  to  modify  the  assigned  errors.  Such 
smoothing  might  consist  of  low-pass  filtering,  spline-function  analysis, 


r 
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piecewise-linear  approximations,  or  global  modeling.  The  present 
analysis  provides  a strong  basis  for  developing  error  assignments  for 
smoothed  data. 

BACKGROUND 

In  this  discussion,  t represents  the  elapsed  time  and  h the  initial 

or  Lagrange  gage  position.  Also,  r(h,t),  v(h,t),  s(h,t),  and  e(h,t) 

represent,  respectively,  the  position,  velocity,  stress,  and  strain  at 

the  point  (h,t).  The  velocity  and  stress  gages  measure  v(h,t)  and 

s(h,t)  at  several  locations,  h1th_,  ....  These  measurements  suffer 

i i 

errors  with  standard  deviations  <jv  and  as,  respectively.  In  addition, 
survey  errors  fih^,  6h2,  6h3,  ...  exist;  these  errors  give  the  distance 
between  the  actual  gage  position  h^  and  the  assumed  one  hm,  m = 1,2,  ... 
M. 

Stress  and  strain  profiles  are  obtained  by  looking  at  the  spatial 
and  temporal  derivatives  and  integrals  of  the  shock-wave  records.  The 
velocity  and  stress  records  are  denoted,  respectively,  by 

v(hm,tn),  s(hm,tn) 


where  m = 1,  2 ...  M,  and  n = 1,  2 ...  N.  The  derivatives  and  integrals 


are  obtained  with  finite-difference  techniques.  For  example,  particle 
position  rOi^t  ) is  calculated  from  the  velocity  samples  using: 


where  6t  is  the  sampling  interval  and  h,,,  is  the  initial  particle  posi- 
tion. Similarly,  the  particle  acceleration  (at  constant  h)  is  obtained 
with: 


— (h  .t  } = v(hm*tn+l)  " v(hm-tn) 

at  m’  rt 


Finally,  the  velocity  gradient  (at  constant  time)  is  approximated  by 
Sv,,  , „ VOlm+l.tn)  - v(hn,,  t„) 


2l(h  t ) _ v(hnd-l»tn)  - vO^.t^ 

3"  n " ^ 


where  Ah  is  the  gage  separation,  here  assumed  constant. 

The  finite-difference  procedure  makes  important  use  of  the  identity: 


<5v  _ dv  £v  At 

dh  t “ dh  path  " ^ h Ah  path 


When  errors  exist  in  the  right-hand  quantities,  it  is  important  to 
select  a path  line  that  minimizes  the  error  assigned  to  the  left-hand 
quantity.  As  discussed  by  Grady7  and  Seaman,18  the  path  line  that  con- 
nects corresponding  features  of  the  gage  records  tends  to  minimize  such 
errors.  Accordingly,  this  path  is  used  in  what  follows  (see  Figure  D-l). 

For  simplicity,  the  following  (unambiguous)  notation  is  used: 


_ hi 

hi 

at  h at 

ah 

dvl 

dv 

At 

dh 

) path 

dh 

Ah 
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FIGURE  D-1  CONSTANT-FEATURE  PATH 
Note  that  tan  (CH  At/Ah. 
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To  avoid  excessively  cumbersome  expressions,  all  errors  will  be 


presumed  uncorrelated,  even  when  it  is  clear  that  they  are  correlated. 
This  procedure  may  decrease  the  error  assignments  as  much  as  50%,  or  it 
may  significantly  increase  them.  In  either  case,  a reasonable  upper 
bound  is  obtained. 

Finally,  the  errors  introduced  by  the  finite-difference  analysis 
are  ignored.  Such  errors  can  be  quite  significant,  especially  on  partial 
derivatives  across  two  gages,  e.g.,  dv/^h. 

STRAINS 

The  tangential,  radial,  and  volumetric  strains,  are  defined  by 

e t = 1 - (r/h) 

ep  = 1 - (dr/dh) 

2 

e = 1 - (r/h)~  (dr/dh) 

V 

where  r is  given  above.  These  quantities  are  estimated  from  the  veloc- 
ity data  using  finite-difference  techniques.  The  errors  introduced  in 
er,  e^.,  and  ey  by  amplifier  and  position  errors  are  considered  in  turn 
below. 

Tangential  Strain 

The  tangential  strain,  e^,  depends  linearly  on  the  position  r. 
Accordingly,  only  this  quantity  is  considered  in  what  follows. 

Amplifier  Noise.  Let  the  amplifier  noise  at  the  gage  position  h 
at  the  sample  time  t be  denoted  by  q(h,t).  Then  the  position  estimate 
will  have  the  basic  form: 

r(*1m,tn^  = *'m  + ^ [v(*1m>'*-i^  + fit 


D-5 


The  statistical  variation  in  the  position  estimate  is  then  given  by: 


! 


l 

1 


k 


V[r]  = E[ (r-r)2] 

= e[{£  qCh.ti)  6t}2] 

2 , .2 

= n ctv  6 t 

2 

where  qv  is  the  variance  of  the  noise.  The  variance  in  the  tangential 
strain  is  thus 


V[r/h]  = nay2  6t2/h2 

Because  the  tangential  strain  involves  an  integral  over  the  experi- 
mental data,  amplifier  noise  has  only  nominal  impact  on  the  estimated 
value  for  the  tangential  strain. 

Position  Errors.  Let  6h  equal  the  difference  between  the  true  gage 
location  h'  and  the  surveyed  one  h.  Then  the  r estimate  will  suffer 
error  because  the  velocity  data  taken  at  the  position  h'  will  be  used 
to  estimate  r at  h,  where  h'  = h + gh.  That  is, 

r(h,ta)  = h + £ v(h' .tj)  6t 

i 

= h + £ v(h  + 6h,t1)  fit 
i 

= h + £ vGi.tp  fit  + £ (-h-,t*)  gh  St 

i i Sh 


I 
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The  first  two  terms  represent  the  usual  finite-difference  formula  for 
the  position  r.  The  last  term  is  an  error  term.  The  variance  of  this 
term  is: 


V[r]  = E[  (£  SX  6h6t)2] 
i 9h 

To  evaluate  this  expression,  the  summation  must  be  writtern  in  terms  of 
the  measured  data. 

The  above  error  assignment  contains  gradients  with  respect  to  h. 
These  gradients  are  evaluated  along  a constant-feature  path  using  the 
identity: 


9v  _ dv  Sv  it 
dh  dh  ?)t  Ah 

The  ratio  At/Ah  can  be  extracted  directly  from  the  experimental  data. 
(This  extraction  gives  At,  since  Ah  comes  from  the  experimental  layout.) 

The  particle  acceleration  is  obtained  by  differentiating  the  veloc- 
ity record: 


dv 

at 


v(hm-tn+l)  - v(hm’tn) 
6 t 


Similarly,  the  total  derivative  is  obtained  by  differentiating  across 
two  velocity  records: 

*1  = v<Vi.  tn+  At>  - V<W 

dh 

Ah 

The  last  three  expressions  allow  evaluation  of  the  variance  in  the 
r-estimate  due  to  position  errors. 
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Because  particle  velocity  varies  slowly  with  h,  gage  position 


errors  will  have  only  nominal  impact  on  the  tangential  strain. 

Radial  Strain 

The  radial  strain,  er,  depends  linearly  on  the  partial  derivative 
gr/^h.  Accordingly,  only  this  quantity  is  considered  in  what  follows: 

Amplifier  Noise.  The  position  gradient  is  obtained  from  the 
constant-feature  path  using  the  identity: 

ar  = dr  ar  At 
?)h  dh  $t  Ah 

The  particle  velocity  Qr/^t  at  constant  h is  an  experimental  measurement 
v(hm,tn),  with  standard  error  The  total  derivative,  on  the  other 

hand,  depends  on  two  velocity  re:ords: 


= r(hm+l>  V At)  - r(h‘m>tn) 


the  standard  error  in  the  r-estimate  was  obtained  above  in  +he  discussion 
of  amplifier  noise  for  tangential  strain. 

The  gage  records  are,  of  course,  statistically  independent.  Also, 
we  are  ignoring  incidental  correlations  among  quantities  derived  from 
the  same  data  set,  e.g.,  rO^,^)  and  v(hm,tn).  The  variance  in  gr/gh 
then  takes  the  uncorrelated  form: 

V[^r/^h]  = 2v[r]/Ah2  + (nv  At/Ah)2 

(The  V[r]  was  obtained  above.)  Generally  speaking,  the  second  term 
dominates  the  first.  Since  At/Ab  « 1,  the  noise  in  Qr/gh  is  comparable 
to  that  in  the  v measurements. 
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Position  Errors.  The  position  error  will  cause  a corresponding 
error  in  the  partial  derivative  Br/Bh.  In  particular, 


if  = V0,,’tn)  = V(U'tn>  5h 


The  partial  derivative  Bv/Bh  was  estimated  in  the  discussion  of  tangential  strains 


above . 


Position  errors  6h  also  propagate  to  the  total  derivative. 


IE  = r<h,m+l.  tn+At)  - r(h  tn) 
dh  


Ah 


_ r<hm+l»t„+At>  “ r(h  ,t  ) 

= — o ju a_ 


Ah 


a^hm  + l*  tn+At>  6lWl 


3h 


Ah 


^r(hm'tn)  6lV 


9h  Ah 


The  last  two  terms  represent  the  first-order  impact  of  position  errors 
on  the  total  derivative.  These  errors  are  statistically  independent, 
since  6hm  + ^ and  6*^  are  independent. 

The  variance  in  Br/Bh  due  to  position  errors  thus  takes  the  form: 


H «Jl\2 

ah]  \Bh  Ah)  ' ah  Ah  5 / 


In  practice,  the  second  term  dominates  the  first  since  r is  an  integrated 
quantity,  while  v is  not. 


Volumetric  Strain 

The  volumetric  strain,  ev>  depends  quadratically  on  the  position 
r and  linearly  on  the  partial  derivative  Br/Bh.  Thus,  the  total  error 
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consists  of  two  terms: 


fi)  =fa)1 


dr  + 2r  dr 
dh  h a*1 


Mr) 


Compared  with  ^r/^h,  the  errors  in  r are  quite  small.  Thus,  we  con- 
centrate on  the  first  term  in  what  follows. 

Amplifier  Noise.  The  amplifier  noise  affects  the  dr/dh  estimate 
and  thus  also  the  volumetric  strain.  The  variance  in  the  volumetric 
strain  is,  thus,  approximately 


r r 2 csrl 

r2  far] 

bid 

h2  l ahj 

where  V[^r/gh]  is  given  above  in  the  discussion  of  radial  strain. 

Position  Errors.  Position  errors  affect  the  Sv/^h  estimate  and 
hence  also  the  volumetric  strain.  In  particular, 


fr2  Sri  r2  I dr  I 

vb  id = i?  vlid 


The  variance  V[av/ah]  was  derived  above  in  the  discussion  of  radial 
strain. 

STRESS  DIFFERENCE 

The  stress  difference  $ depends  both  on  the  velocity  and  stress 
records.  In  addition,  this  quantity  is  obtained  as  a small  difference 
between  two  large  quantities.  That  is,  $ has  the  approximate  form: 


2 l at  ah/ 

(The  two  partial  derivatives  have  opposite  signs.)  Typically,  the 
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magnitude  of  $ is  10%  of  either  dv/^t  or  Qs/^h.  Accordingly,  small 
errors  in  either  partial  derivative  lead  to  large  errors  in  the  stress 
difference.  The  situation  is  further  aggravated  by  the  fact  that  both 
partials  are  obtained  as  derivatives  of  the  experimental  data,  a cal- 
culational  process  that  also  magnifies  noise. 

Amplifier  Noise 

The  amplifier  noise  in  the  velocity  gages  loads  to  a noisy  accelera- 
tion estimate.  That  is, 


9v  _ v(h,tn+1)  _ V(h,tn) 
9t  - 

lias  a variance  of  the  form: 


For  small  6t,  this  variance  becomes  quite  large.  Conversely,  large  §t 
reduces  the  variance,  but  also  increases  the  finite-difference  errors. 
The  stress  gradiant  is  obtained  with  the  path-line  technique: 

5s  ds  9s  At 
dh  dh  J)t  Ah 

Tlie  variance  in  the  first  term, 


is  just 


ds 

dh 


s%+i,tn+^  - styy 

Ah 
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Similarly,  the  second-term  variance  is 


The  total  variance  due  to  amplifier  noise  follows  by  addition: 


We  are  again  ignoring  minor  correlations  that  occur,  for  example,  be- 
tween ds/dh  and  3s/9t.  In  practice,  At  » Ah  » fit,  so  that  the  first 
and  third  terms  dominate  the  second. 


Pos it ion  Errors 

Position  errors  fih  affect  the  partial  derivatives  Qs/gh  used  to 
calculate  the  stress  difference  §.  The  particle  acceleration  5v/dt  will 
be  in  error  according  to 

*1  = v<hVWl)  ~ v<»  B'V 

_ 

= v(VW  - vOVi.tn) 

6 1 

+ ?v(hm>tn  + l ) i hm  _ c>v  (hm>  tn) 

Bb  6 1 ^h  6 t 

Tlie  last  two  terms  represent  the  induced  errors;  they  can  be  written 
2 

as  (a  v/ahat)6h.  The  variance  in  the  acceleration  estimate  is  thus 


2 


The  second  derivative  can  be  estimated  with  the  formula  given  in  Appendix  E. 
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The  position  errors  6h  also  affect  the  stress  gradient: 


9s  _ ds  5s  At 
9h  dh  £t  ^h 

By  symmetry,  the  variance  in  9s/9t  is  given  by: 


The  second  partial  is  obtained  by  the  formula: 

52v  _ [Sv] 

5t9h  9h  L^tJ 

= £_  r^i  _ a_  rizi 

dh  LstJ  " 9t  LstJ  Ah 

= <v(hm+l,tn+i+At)  - vdim+t.tn+At) 

' 6t  Ah 


“ v(hm,tn+l)  - vOytn)  ) 

6t  At  ) 

- jvCWi+M)  ~ ^W^t) 

l 6t  At 

v(hm>tn+l)  - vCh^jtjj)  \ 

6t  At  J 

The  total  derivative  has  the  form: 

= s(h,m+i*tn+At)  - s(h',tn) 
dh  ~ 

= s(hm+1,tn+At)  - s(hm,tn) 

Ah 

9s  , Ss  , 

+ & 6h™+1"^  6h™ 
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Here  the  two  6h  are  uncorrelated.  The  total  variance  in  ds/dh  is  thus 


where  Qs/^h  is  calculated  in  the  usual  way. 


The  total  variance  in  the  stress  difference  is  thus 

In  many  instances,  the  second  derivative  will  prove  mild,  thus  allowing 
the  second  term  to  dominate  the  error  assignment. 


Appendix  E 

STRESS-STRAIN  ERROR  ASSIGNMENTS 
FOR  GLOBAL  LAGRANGE  ANALYSIS 

D.  G.  Falconer 


INTRODUCTION 

Global  Lagrange  analysis  of  shock-wave  field  data  offers  several 
advantages  over  finite-difference  techniques.  In  particular,  global 
analyses : 

• Reduce  the  impact  of  random  amplifier  noise  and 
gage  position  errors  by  systematically  averaging 
(smoothing)  the  experimental  data. 

• Allow  enforcement  of  known  physical  constraints 
(e.g.,  causality)  during  the  fitting  effort. 

The  chief  purpose  of  this  appendix  is  to  develop  error  assignments  for 

the  stress  and  strain  trajectories  calculated  by  the  global  Lagrange 

analysis . 


BACKGROUND 


Shock-wave  experiments  generate  velocity  v(h  ,t  ) and  stress 

m n 

s(h  ,t  ) data  at  various  gage  locations  h ,m  = 1 , 2,  3 . , . M.  Each  of 
m n m 

these  may  be  digitized  at  regular  times  t , n = 1 , 2 , 3 . . . N.  In  a 

n 

global  fitting  procedure,  these  space-time  profiles  are  used  to  determine 
a best  set  of  model  parameters  a , j = 1 , 2 , 3 . . . R . In  particular,  the 
fitting  procedure  attempts  to  minimize  the  chi-square  error  between  the 
experimental  values  and  model  function  f, 


f = f(h,t,k;  a^ ) 

Here  k = 1 for  velocity  data,  and  k = 2 for  stress  data. 
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The  chi-square  error  has  the  basic  form: 


V2  = Z Tf(h  ,t  ,l;a  ) - v(h  ,t  )1 2/a2 
m n j m n v 


+ E [f  (h  ,t  ,2;  a ) - s(h  ,t  )]2/(j2 
m n j m n s 


In  practice,  the  standard  errors  in  v and  s may  also  depend  on 
position  and  time,  that  is, 


a = O (h  ,t  ) 
v v m n 


a = a (h  ,t  ) 

s s m n 


However,  adjacent  measurements,  i.e.,  at  (h  ,t  ) and  (h  , t ) are 

m n m n+1 

presumed  uncorrelated. 

The  chi-square  fitting  algorithm  minimizes  chi-square  by  picking 

a best  set  of  numerical  values  for  the  parameters  a . During  the 

j 

fitting  effort,  the  computer  algorithm  also  obtains  an  error  covariance 

matrix  for  the  parameters  a . This  matrix 

j 


Via J = E[ 5a  5a  ] = [A  ], 

i j i j 


gives  the  variances  (A  ) and  correlations  (A  ) in  the  fitted  values 

11  i j 

a^.  (Note  that  A is  a single  matrix  element,  not  an  algebraic  product.) 

Once  a best  fit  is  obtained,  the  function  f is  used  to  deduce 
stress  and  strain  profiles.  In  particular,  the  following  quantities 
are  of  interest: 
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e = 1 - r/h 

e = 1 - dr/dh 
r 

e = 1 - (r/h)2  (dr/3h) 
v 


tangential  strain 
radial  strain 
volume  strain 


stress  difference 


t 

r(h,t)  = h + f v(h,t/  )dt 7 particle  position 

o 


(The  derivative  notation  used  above  is  explained  in  Appendix  D.) 
Accordingly,  the  principal  objective  of  this  appendix  is  to  use  the 
errors  assigned  to  the  parameters  a . to  develop  errors  for  the  stress 


and  strain  quantities. 


J 


STRAINS 


For  simplicity,  we  write 


f = f (h,  t , 1 ; a ) 

v 0 


f = f (h,t,2;a  ) 
s J 


Strain  estimates,  of  course,  make  use  of  only  f . 


For  definiteness,  we  assume  the  following  forms  for  f and  f : 

v s 


f =(ax+ax  + a x )(i  + a h + a h ) 
v 1 2 3 4 5 


2 3 2 

f =(ax+ax  +ax)(l+ah+ah) 
s 6 7 8 4 5 


Here  x represents  the  elapsed  time  since  shock  arrival  at  the  gage 
position  h: 


x = t + a h 
9 

where  aq  is  always  a negative  number.  This  nine-parameter  model  meets 

causality  if  f = f =0  for  x < 0. 
v s 
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The  particle  position  r(h,t)  is  used  for  both  the  tangential  and 


radial  strains.  It  has  the  explicit  form: 


r(h,t)  = h + f f (h.t^dt 

V 


2 3 4 

x x x 2 

= h + (a  — + a — +a  — )(1  + a h + a h ) 
12  23  44  4 5 


Small  variations  6a  in  the  fitted  parameters  produce  a small  variation 


6r  in  the  position  estimate: 


6r=r|i  6a. 
S£,j  J 


= *r  6ai  + T 6a2  * T 6*3H1  * V + a5h  > 

+ "ir1!  r’ <h  5a4 * "X* 


The  variance  in  the  position  estimate  is  thus  given  by; 


vm  - * r fr A . 

. da  da  ij 
i , J i J 


where  the  quantity  A comes  from  the  fitting  program. 

ij 

Tangential  Strain 

The  tangential  strain,  et»  depends  linearly  on  the  position  r. 
Thus,  the  assigned  error  follows  readily: 


v[e  ] = V[r]/h 
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Radial  Strain 


The  radial  strain,  er , depends  linearly  on  the  partial  derivative 
dr/dh.  According  to  our  model: 


— _ 1 , x X x 

ah  = + (ai  r + a2  r + a3  (a4  + 2a5h> 


+ (a  x + a x 


2 3 


1 + a2X  + V > a9(1  + a4h  + agh  ) 


Small  changes  fia^  thus  produce  the  following  change  in 
derivative : 


the  partial 


/dr \ 92r 

6 V3h)  = E 9a  dh  6a  i 


J J 


The  variance  in  Sr/dh  and  hence  also  e thus  has  the  form: 


vre  1 = £ -A2r  . 

r i,j  da. Sr  da.dh  ij 


where  A_  obtains  during  data  fitting. 


Volumetric  Strain 


Th.  volumetric  strain,  depends  ou.dratlcally  on  r and  linearly 
on  dr/ah.  Small  change.  in  these  quantities  lead  to  the  following  change. 


6e  = - £-  6 _ 2r  dr 

v h2  W h2  ah 
2 

- - — y d r 2r  dr  dr 

h2  SajSh  aj  ~ h2  sh  E aT  6aj 

= - £ fr  -?2r  dr  drj 

h2  J [ sajdh  dh  3aJ  6aj 
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where  we  have  propagated  through  to  the  parameter  errors.  The  variance 
in  volumetric  strain  is  thus: 

2 

vrf  1 = ^ [ ...  1.  [ ...  1 a. . 

v h4  ij  i j xj 

where  bracket  is  written  out  just  above.  The  respective  partial 
derivatives  are  specified  above  and  in  the  discussion  of  radial  strain. 

STRESS 

The  stress  difference  <f>  depends  linearly  on  dv/Bt  and  5s/Sh.  The 
particle  acceleration  depends  on  the  fitted  parameters  in  the  following 
way 

A 3f 

dv  v 

at  = at- 

Small  changes  in  the  fitted  parameters  thus  lead  to: 


(I) 


a2f 

— — 5a 

Sa  .at  j 


Similarly,  the  stress  gradient  has  the  following  form: 


* Bf 

ds  s 


dh  Sh 


Small  changes  in  a thus  produce  the  following  changes  in  the  stress 
j 


gradient : 


(I)  ' 


a2f 


- — 3a 

Sa^dh  j 


The  above  expression  is  messy,  but  easy  for  a digital  computer. 
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Small  changes  in  the  fitted  parameters  thus  impact  the  stress 
difference  in  the  following  way: 


6 


(I)  * 6(l)) 


. .2  2 

h ^ d v + o s 

2 da  dt  da  dh  j 
j J J 


The  variance  in  the  stress  difference  thus  takes  the  basic  form: 


= (!)%.  c ...  1.  [..•].  A 

\ 2/  i,J  l J lj 


OTHER  MODELS 


Other  models  for  the  experimental  data  are  handled  in  a completely 
analogous  way.  One  such  model  is  the  following: 


2 3 5 

f = Ta  x + a x + a x "|  exp(a  x)/(h+l) 
v i 2 3 k 4 


2 3 5 

f = [ax  + a x + a x 1 exp  (a  x)/(h+l) 
s 6 7 8 4 


where 


x = t + ah+a  h 
9 10 


The  above  formulas  allow  for  exponential  release,  power-law  amplitude 
loss,  and  nonuniform  propagation  speeds. 

Note  that  global  fitting  avoids  the  need  for  path-line  derivatives 
to  evaluate  the  gradients  dr/dh  and  ds/dh. 
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Appendix  F 

STRESS-STRAIN  ERROR  ASSIGNMENTS 
FOR  PIECEWISE- LI NEAR  LAGRANGE  ANALYSIS 

D.  G.  Falconer 

INTRODUCTION 

Previously  in  this  program  we  have  extracted  stress  and  strain 
trajectories  from  field  data  using  the  piecewise-linear  Lagrange 
analysis  (Appendix  A).  This  analysis  consists  of  the  following  steps: 

(1)  The  key  features  of  each  velocity  and  stress  history 
are  identified,  located,  and  measured.  (Both  time  and 
magnitude  errors  exist  with  this  procedure.) 

(2)  Corresponding  features  on  the  several  velocity  and  1 

several  stress  records  are  then  fit  with  exponential- 

type  functions  F^(h)  and  G (h) , respectively.  ^0 

gives  the  particle  velocity  associated  with  the  n feature 

at  the  gage  position  h.  and  G (h)  gives  the  stress  level 

th  n 

associated  with  the  n feature  at  the  gage  position  h. 

(3)  S^jnilarly  functions  T (h)  are  fitted  to  the  times  at  which  the 

n feature  occurs  at  the  gage  position  h.  The  functions 

F , G , T provide  continuous  estimates  of  particle  velocity, 
n n n th 

stress  level,  and  occurrence  time  of  the  n feature. 

Velocities  and  stresses  between  the  key  features  are 

interpolated  linearly  using  F and  G . 

n n 

The  principal  objective  of  this  appendix  is  to  develop  error  assignments 
for  the  strains  and  stresses  derived  from  the  experimental  data  using 
the  piecewise-linear  Lagrange  analysis.  Mathematical  approximations 
are  used  liberally  and  correlated  errors  are  properly  accounted  for. 


BACKGROUND 

Let  the  gage  positions  be  denoted  by  h ,m  = 1,  2,  ...  M,  and  the 
th  m 

occurrence  time  of  the  n feature  be  denoted  by  t . Thus,  following 

n 

step  (1)  above,  the  velocity  and  stress  data  are  of  the  form: 


v(h  , t ) , s (h  , t ) 
m n m n 
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The  standard  errors  in  the  velocity  and  stress  measurements  are  denoted 


by  (j  and  O'  , respectively.  Also,  the  standard  errors  in  the  gage 
v s 

position  and  occurrence  time  are  denoted  by  o and  <7  , respectively. 

h t 

The  fitting  step  (2)  defines  the  parameterized  functions 

, n , n n n 

F (h,  a ) and  G (h,b  ),  where  a and  b are  the  fitted  parameters  for 
n ..i  n i i i 

t h 

the  n features.  During  fitting,  error  covariance  matrices  for  the 

a and  b are  developed.  For  convenience,  these  matrices  are  written 
i i 

in  the  following  form: 


r n _ r n „ n r n 
VTa  1 = ET6a  6a  1 = [A  ] 
i j i j 

V[b  ] = E[6b  6b  ] = TB  1 
i J ij 

n n 

(Here  A and  B are  matrix  elements  (numbers),  not  algebraic  products.) 
ij  i j 

The  diagonal  elements  represent  the  variances  in  the  fitted  parameters; 
the  off-diagonal  ones  represent  the  convariance  between  the  fitted 
parameters . 

r n.  r n-i 

The  error  matrices  V[a  J and  VI b J depend  on  the  standard  deviations 

n , <7  , and  a . The  error  represented  by  rr  is  an  amplitude  error  and  thus 
v s h v 

handled  in  the  usual  way  during  fitting.  The  error  represented  by  q , 

h 

however,  must  be  propagated  to  look  like  an  amplitude  error.  An  easy 

way  to  accomplish  this  task  is  to  first  fit  F (t,a  ) to  the  data  set 

n i 

assuming  =0,  and  then  use  the  fitting  function  in  the  following 


propagation  formula: 


The  errors  represented  by  g and  are,  of  course  incoherent,  so 
that  their  variances  add: 


The  fitting  process  is  then  repeated,  using  a in  place  of  o' 

v 

A similar  approach  is  taken  with  fitting  G (h,b  ).  In  practice, 

n i 

O . CT  , and  o may  depend  on  the  gage  position  h and  the  feature 
v s h m 

number  n. 


Under  step  (3),  the  function  T (h,c  ) is  fit  to  the  occurrence 
. n i 

t h 

times  of  the  n feature.  This  fit  suffers  uncertainty  (oO  in  the 

time  estimate  and  uncertainty  (o  ) in  the  position  estimate.  The 

h 

first  uncertainty,  being  an  amplitude  error,  is  handled  in  the  usual 

way  during  fitting.  The  second  uncertainty  is  accounted  for  following 

an  initial  fit  with  o,  = 0.  The  fitted  function  is  then  used  to  derive 

h 

a total  amplitude  error: 


In  practice,  both  o and  o may  depend  on  gage  position  h and  the 
t h m 

feature  number  n.  In  other  words: 


CTt  = CTt(m,n) 
ffh  = (Th<m,n) 

This  dependence  is  readily  handled  by  chi-square  fitting  routines. 
STRAINS 

The  tangential,  radial,  and  volumeric  strains  have  the  following 
forms : 

et  = 1 - r/h 

e = 1 — dr/dh 
r 

e = 1 - (r/h)2  (dr/5h) 
v 

F-3 


The  error  assignments  for  these  quantities  are  developed  below. 


Tangential  Strain 

The  main  ingredient  of  the  tangential  strain,  , is  the  particle 
position  r.  This  quantity  obtains  in  the  usual  way  by  integrating 
the  smoothed  (i.e.,  fitted)  velocity  records: 


r(h,t)  = h + f vih.t'ldt7 


= h + — Z Tl  (h,a  ) - F (h,  a )l(t  - t ) 
2 n i n-1  i n n-1 

n=l 


= h + E F (h,a  ) £ t 
n i n 

n 


We  have  introduced  the  notation 


At  =1/3  (2t  - t - t ) 

n n n+1  n-1 


The  occurrence  times,  t , are  taken  from  the  fitted  values,  that  is, 

n 

t (h)  = T (h,c  ).  The  final  time,  t,  is  interpolated  if  it  falls 
n n i 

between  two  key  features. 

n 

Small  changes  6a  in  the  fitted  parameters  then  lead  to  the 
following  changes  in  the  particle  position: 


_ n 2 n 

Ar  = Z A t 6a 

. , n n i 
n,  i da 

i 


(For  simplicity,  we  have  presumed  that  temporal  errors  are  small.) 
The  variance  in  the  r-estimate  follows  immediately 


.*»  v -*  • • - 


Jfl 


IdF  3f 

S -2-  A2t  A2t  6an  6am 

„ m n n m m i j 

M S.(  8. 


SF  9f  2 

- S — — (A2t  ) A". 

n,i,j  3a”  3a”  n lj 

i J 

(The  variations  fia^  and  6a^  are  presumed  statistically  independent 

unless  n = m.)  The  above  formula  is  readily  handled  with  a digital 
computer. 

Radial  Strain 

The  main  constituent  of  the  radial  strain,  e , is  the  partial 

r 

derivative  Sr/3h.  This  quantity  obtains  from  the  r-estimate: 


Sr  n 2 

-s.  = 1 + £ TT-  A t 

oh  n 

n 

Small  variations  in  dF^dh  thus  lead  to  the  following  change  in  the 
partial  derivative: 


For  simplicity  we  presumed  that  the  various  AOF^/Sh)  are  statistically 
independent.  In  this  case  the  variance  follows  directly: 


-*„1  , fSF  9F 

V ~ = E (42t  ) E -2  -s 

on  n,m  n 3h  3h 


Thus,  we  need  the  expectation  of  the  various  products  Of  /3h)  (3f  /3h) 

n m 


The  derivative  dF  /dh  obtains  through  the  path-line  formula: 
n 


dF  dF  dF 
n n n dt 


dh  dh  dt  dh 


The  ratio  dt/dh  follows  directly  from  the  fitted  times  T (h,c  ): 

n i 


dT  (h,c  ) 
dt  n i 

dh  - dh 


= tan  (lb  ) 
n 


Similarly,  the  total  derivative  dF  /dh  follows  immediately  from  the 

n 

fitted  velocity  F . The  time  derivative  is  obtained  by  interpolating 
n 

between  key  features. 


n+ 1 n 

W y1|h-,1  1 ~ V'V 

St  = T ,(h,cn+l>  - T <h,c"> 
n+1  l n i 


In  what  follows  the  denominator  is  denoted  by  At  . 

n 

We  presume  that  the  principal  source  of  strain  error  results  from 


F , rather  than  from  T . In  this  case  small  changes  in  the  parameters 
n n 

an  lead  to  the  following  changes  in  the  partial  derivative: 
i 


= z 

\dhy 

i 

aFn+l 

. n+1 
6s, 

dF 

n 

da°+1 

At 

n 

da" 

6an\ 

- — — - It  an  (ib  ) 

n At  I n 


F-6 


The  expectation  of  the  cross  partials  is  thus: 


SF 

dF 

n 

m 

dh 

dh 

= E 


E jrn.il  [m, jl  j 
i.j  ' 


where  the  bracketed  quantities  abbreviate  the  bracketed  quantities  in 
the  penultimate  expression. 

. n _ m 

The  variations  oa  and  6a  are  statistically  independent  unless 

t j 

m = n.  If  n = m,  the  last  expression  becomes: 


dF 

n 

dF  ’ 
n 

- r( 

±11 

F \ dF 
nl  n 

tan  lb 

n 

dh 

dh 

L i 

i J n 

' da 

i 

At 

n 

_ 

_a.| 

<d\) 

dF  tan  tb 
j n n 

,»  , 8!»«  ®Vl  (*“  V»\.->  > 

n * 

da 

L J 

idh  ) 

' * n At 

da  n 

J -1 

11  a.""1  a."*1  ' > ,J  1 

i j 

when  n = m + 1,  then  we  have  instead: 


dF 

n 

dh 

dF  1 

n-1 

dh 

= E 

'jL^y 

_ n l dh 

dF 

^ n 

tan  ib 

n 

At 

- 

i . J 

da  ' 

L i 

da 

i 

n 

dF  tan  >/, 

n 1 n-1 

. n At  , 

da  n-1 
L j 


i j 


If  n + 1 = m,  we  get  a similar  expression  with  n everywhere  replaced 
by  n + 1 : 


dF 

n+1 

dF 

n 

= E 
i.j 

d 

dF  tan  lb 

n+1  n+1 

dh 

dh 

, n+1 
da 

i 

It) 

(dh  j 

n+1  At  , 
da  n+1 

i 

dF  , tan  I b 
n+1  n 


n+1  At 
da  n 


vn+l 

U 


All  other  expectations  are  zero. 


Volumetric  Strain 


The  volumetric  strain,  e , depends  quadratically  on  r and  linearly 
on  Or/dh)  . Accordingly,  small  changes  in  these  quantities  produce  the 
following  changes  in  the  volumetric  strain: 

- (iMs)  ♦ (?)  (i)  - 

Since  r represents  an  integrated  quantity,  the  error  in  dr/dh  is 

expected  to  dominate  the  error  assignment  f or e . In  other  words, 

v 

the  volumetric  strain  has  the  following  variance: 


The  right-side  variance  was  obtained  above. 

STRESS  DIFFERENCES 

The  stress  difference,  $,  has  the  approximate  form: 


_ h Bv  ds- 

* ' 


The  estimation  procedures  for  the  two  derivatives  are  statistically 
independent.  Accordingly,  the  variance  in  $ has  the  form: 

j bl’ * ''tfi. 


Thus,  we  need  the  two  variances  on  the  right-side  expression. 
The  time  derivative  is  approximated  in  the  usual  way: 


* 

i 


. 


where  n is  selected  so  that  T < t < T . Small  changes  in  the 

n n+1 

parameter  values  then  lead  to  the  following  changes  in  the  time 
derivative: 


At 

n 


'Thus  the  variance  has  the  form: 


1 

/ 3F 

SF  \ 

3v 

I n+1 

n+1  \ 

St 

= £ 

l ^ n+1 

^ n+l  / 

i.  j 

' Sa 

Sa  ' 

i 

j 

. n+1  n 

since  oa . and  oa  are  statistically  independent. 


The  space  derivative  is  estimated  in  the  usual  way  as  well: 


- dG 

3s  ds  os  dt  n 

Sh  dh  St  dh  dh 


G - G 
n+1 n 

At 


tan  ij) 


Small  changes  in  the  fitted  values  then  lead  to  the  following  change  in 
the  space  derivative: 


The  variance  of  Ss/Bh  then  follows  immediately: 


In  practice,  the  indicated  derivatives  are  evaluated  numerically  and 
the  sum  handled  as  matrix  multiplications. 
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